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FOREWORD 


This  report  presents  results  of  a study  conducted  by  Westinghouse 
Electric  Corporation,  Baltimore,  Maryland,  under  Air  Force  Contract  F33615 
77-C-2074,  "Integrated  Thermal  Avionics  Design"  during  the  period  of 
September  1,  1977  to  June  1,  1978.  This  contract  was  conducted  under  the 
sponsorship  of  the  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson 
Air  Force  Base,  Ohio  with  Dr.  A.  H.  Mayer  (FEE)  as  project  engineer. 


The  results  of  the  study  and  recommendations  for  further  study  are 
presented  herein.  The  Westinghouse  Program  Manager  was  R.  F.  Porter. 
Other  major  contributors  were  R.  T.  Dolbeare,  E.  R.  Levitt,  R.  H.  Warsham 
R.  J.  Bennett,  J.  T.  Prichard,  Y.  Lord,  R.  A.  Marlett. 
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INTRODUCTION 


The  primary  goal  of  Integrated  Thermal  Avionics  Design  (ITAD)  is  the 
development  of  a system  of  software  programs  to  accomplish  optimum  life 
cycle  cost  (LCC)  design  of  avionics  systems.  Accomplishing  this  goal  means 
ultimately  making  all  necessary  trade-off  decisions,  not  only  concerning  the 
avionics  system  design  itself,  hut  also  concerning  the  aircraft  environ- 
mental control  system  (ECS). 

This  final  report  details  the  studies  and  progress  made  toward  these  goals 
within  the  constraints  of  this  initial  phase  of  the  project.  This  initial 
phase  was  not  intended  to  yield  a working  software  system.  It  is  a study 
phase  intended  to  design,  in  an  architectural  sense  only,  the  software 
system  required  to  accomplish  LCC  optimization. 

One  of  the  major  ground  rules  of  this  study  is  to  use  existing  software 
wherever  possible.  Thus  a survey  of  existing  computer  programs  has  been 
made  to  determine  which  are  most  applicable  to  the  ITAD  system.  In  the 
initial  ph„se  of  the  program,  it  was  recommended  that  the  ITAD  software 
should  be  designed  for  eventual  implementat ion  on  the  CDC-6600  computer 
located  at  Vs  right -Patterson  Air  Force  Base  (WTAKBl.  Therefore,  part  of 
the  choice  of  computer  programs  is  based  on  their  applicability  to  this 
hardware  facility  as  well  as  their  applicability  to  the  ITAD  system. 

Many  functions  required  to  be  performed  by  the  LCC  opt imi zat ion  process  or 
to  accomplish  some  peripheral  tasks  are  non-existent.  In  these  cases, 
algorithms  have  been  defined  to  fill  these  gaps.  Software  code  lias  not 
been  written  in  this  phase  of  the  study  except  where  it  was  convenient  or 
made  completion  of  this  study  easier. 

Three  separate  groups  of  effort  are  reported  herein.  These  are  divided 
into  two  major  tasks  with  one  task  having  two  separate  phases.  Task  1 
involves  the  definition  and  design  of  the  ITAD  software  system.  Task  II 
involves  an  actual  reliability  optimization  of  an  existing  pair  of  Printed 
Wiring  Boards  (PWB)  and  predictions  of  the  influence  this  optimization 
process  has  on  total  LCC. 

Although  optimization  of  LCC  for  every  type  of  electronic  component  and 
subassembly  is  desirable  and  eventually  possible,  it  became  immediately 
apparent  that  effort  should  initially  be  concentrated  on  the  digital  PWB 
area.  Since  the  majority  of  parts  in  modern  avionics  systems  are  contained 
in  the  digital  area,  large  gains  in  total  system  reliability  and  Logistics 
Support  Cost  (LSC)  can  be  realized.  Automated  wire  routing  programs  are 
best  suited  for  wiring  digital  circuitry.  Also,  thermal  optimization 
techniques  are  relatively  straight  forward  since  modern  digital  packaging 
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techniques  allow  consideration  of  optimum  arrangement  on  a row  by  row  basis. 
Interactions  between  rows  are  a second  order  effect  for  most  packaging 
techniques  commonly  used.  In  contrast,  most  analog  boards  (or  boards 
consisting  of  discrete  devices  instead  of  microelectronic  devices  (MED's) 
have  very  complex  arrangements  requiring  computations  of  temperatures  and 
resulting  reliability  based  on  complex  interactions  between  all  the  compo- 
nents on  the  PKB.  In  addition,  changing  the  location  of  one  component  may 
require  a complete  rearrangement  of  all  the  components  in  a board  - probably 
not  in  any  systematic  manner. 

The  optimization  algorithms  developed  in  this  study  for  digital  Printed 
Wiring  Assemblies  (PWA)  have  the  necessary  flexibility  that  they  can  be 
used  to  determine  positions  of  PWA's  relative  to  each  other  within  a unit. 

In  addition,  they  can  be  extended  to  the  placement  of  each  unit  making  up 
the  total  avionic  system.  Initially  this  total  optimization  cycle  can  be 
performed  by  the  manufacturer  of  each  avionic  system  going  into  the 
receptor  aircraft.  Finally,  the  same  procedure  can  be  used  for  the  entire 
aircraft  avionic  systems  either  by  the  airframe  manufacturer,  the  customer, 
or  an  avionics  integration  subcontractor. 
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SECTION  II 


TASK  I - ITAD  DEFINITION 


Task  I of  the  ITAD  study  was  divided  into  two  efforts.  The  first,  to 
recommend  the  hardware  system  on  which  to  implement  the  ITAD  software.  The 
second  to  define  the  ITAD  software.  These  separate  goals  are  entitled  Phase 
I and  Phase  II  respectively  and  are  discussed  in  detail  in  sections  2.1.1 
and  2.1.2. 

2.1  PHASE  I,  COMPUTER  FACILITY  DEFINITION 

The  initial  phase  of  the  ITAD  was  aimed  at  choosing  the  best  computer  site 
for  initial  installation  of  the  ITAD  software.  Three  alternatives  were 
explored : 

1)  The  DEC-10  computer  at  KPAKB  Avionics  Lab 

2)  The  C DC-6600  at  WPAKB 

3)  A third  alternative 

Among  those  systems  considered  for  the  third  alternative  were  the 
Westinghouse  Univac  1110,  a commercial  timeshare  system,  and  a combination 
of  both  WPAFB  systems. 

The  commercial  systems  were  eliminated  due  to  expense  and  because  they 
were  CDC  machines  similar  to  those  at  WPAFB  ASD  and  thus  had  no  advantage 
in  machine  type.  The  Westinghouse  Univac  1110  system  was  eliminated  because 
access  to  the  computer  from  terminals  at  WPAFB  could  not  be  provided. 

The  remaining  computer  sites,  i.e.,  DEC-10,  ASD  CDC-6600,  or  a combination 
of  the  two  were  compared  on  the  basis  of  the  following  criteria: 

1)  Conversion/instal lation  costs 

Most  of  the  software  to  be  used  is  located  on  the  CDC-6600.  Those 
which  are  not  are  easily  converted  from  the  DEC-10  or  is  available 
in  a form  compatible  with  the  6600. 

2)  Availability 

This  is  an  assessment  of  the  availability  of  lack  of  the  candidate 
machines  at  the  given  time  a user  tries  to  start  an  ITAD  run. 

3)  Additional  facilities  required 

This  criteria  considered  the  need  for  additional  equipment  such  as 
graphics  terminals,  microprocessors , available  memory,  etc.  The 
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CDC-6600  system  presently  has  all  the  foreseeable  facilities 
required.  The  inclusion  of  STAGING  on  the  DEC-10  would  require 
at  least  one  additional  terminal. 

4 ) Ease  of  Use 

This  considered  the  amount  of  knowledge  and  experience  required 
by  the  user  to  use  ITAD  on  each  of  the  three  systems  under 
consideration.  A combination  system  would  require  a higher 
level  of  knowledge  and  more  difficulty  for  the  user. 

5)  Response  time 

Response  time  is  defined  as  the  time  between  starting  'he 
problem  and  obtaining  an  answer.  This  is  dependent  on  avail- 
ability. Machine  speed,  loading  of  the  machine. 

6)  Software  Maintainability 

Maintainability  expenses  may  vary  between  the  candidate  sites. 
This  includes  costs  of  updating  as  well  as  eliminating  bugs 
found  during  use. 

7)  Difficulty  in  installation  in  another  facility  at  later  date. 

The  ease  of  installation  of  ITAD  on  other  computer  sites  after 
having  been  installed  at  WPAFB  will  depend  on  the  original  site 
chosen . 

A combination  of  machines  is  felt  to  be  undesirable  because  of  difficulties 
of  intermachine  interactions  with  the  Data  Base  Manager  (DBM)  and  because 
of  increased  maintenance  costs.  Also,  since  the  ITAD  software  will  be  used 
by  various  avionics  manufacturers,  it  will  need  to  be  converted  to  their 
computer  systems.  This  will  be  more  difficult  if  the  original  software  is 
split  between  two  machines. 

Thus  the  combination  of  machines  was  eliminated  as  the  final  ITAD  hardware 
facility.  However,  much  of  the  software  required  for  ITAD  is  currently 
resident  on  Kright-Patterson  computers  but  is  split  between  the  Flight 
Dynamics  Lab  facility  (ASD)  and  the  Avionics  Lab  facility. 

There  are  at  the  present  time  some  facility  problems  with  both  of  these 
sites.  The  DEC-10  at  AFAL  is  dedicated  to  real  time  simulation  causing 
frequent  unacceptable  turn-around  times.  The  AFFDL  CDC-6600  permanent  and 
scratch  file  systems  are  overloaded  causing  frequent  periods  of  inaccessi- 
bility. Elimination  of  these  problems  is  a matter  of  speculation,  but  it 
seems  that  the  most  likely  system  to  overcome  these  problems  is  the  AFFDL 
CDC-6600  since  its  primary  mission  is  to  provide  the  entire  base  and  the 
Flight  Dynamics  Lab  in  particular  with  required  computational  facilities. 
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Although  the  objective  of  this  phase  of  study  was  only  to  choose  a computer 
| site,  a parallel  investigation  was  undertaken  to  determine  the  particular 

software  and  their  interactions  required  to  meet  the  second  phase  goals. 

This  was  necessary  to  make  an  intelligent  judgement  of  the  type  of  system 
and  capabilities  required  of  the  site.  Based  on  this  study,  the  block 
diagram  shown  in  figure  1 was  developed. 

Based  on  these  requirements,  a two-step  implementation  was  decided  to  be 
most  advantageous.  The  wire  routing  program,  PCPRA,  is  currently  resident 
on  the  AFAL  DEC-10  system.  In  addition,  the  personnel  with  the  expertise 
to  use  this  program  (which  does  require  some  interactive  manipulations)  is 
also  at  the  AFAL.  It  is  also  felt  that  some  period  of  study  is  required 
before  this  optimization  loop  will  become  completely  automated.  Therefore, 
the  first  step  of  implementation  was  chosen  to  put  all  required  1'l‘A!)  func- 
tions on  the  ASD  AFFDL  CDC-6600  facility  except  the  wire  routing  program. 
This  allows  the  utilization  of  the  expertise  at  AFAL.  Once  the  best  imple- 
mentations of  the  PCPRA  program  in  the  ITAD  system  has  been  defined  and  AFFDL 
personnel  are  familiar  with  its  use,  the  program  can  be  converted  to  the 
CDC-6600  and  integrated  into  the  ITAD  system  to  run  in  an  automatic  optimi- 
zation algorithm. 

Thus  figure  1 will  initially  have  a break  at  the  input  to  and  output  of 
the  wire  routing  program  representing  a change  of  computer  hardware 
facilities.  In  step  two,  this  break  will  be  eliminated  by  installing  PCPRA 
on  the  CDC-6600. 

2.2  PHASF;  II,  ITAD  SOFTWARE  DEFINITION 

Phase  II  of  Task  I is  dedicated  to  the  definition  and  description  of  the 
software  required  and  their  interactions  required  to  accomplish  the  various 
functions  needed  to  optimize  avionics  designs  for  minimum  I.CC,  perform 
various  environmental  analyses,  and  various  reliability  and  statistical 
calculations.  The  exact  nature  of  these  functions  is  described  in  the 
section.  Wherever  possible,  existing  software  is  used.  The  choice  of 
software  is  biased  toward  those  existing  at  WPAFB  or  its  compatibility 
with  the  CDC-6600  computer  at  WPAFB. 

Where  existing  software  has  not  been  found  to  apply  to  a function  required 
by  ITAD  capability  requirements,  algorithms  are  defined. 

2.2.1  LCC  Optimization  Loop 

For  purposes  of  the  discussion  of  the  optimization  loop,  the  example  will  be 
centered  on  the  optimization  of  the  Printed  Wiring  Assembly  (PWA)  type  shown 
in  figure  2 which  will  be  referred  to  as  "coplanar  cooled."  This  assembly 
consists  of  the  boards  and  heat  is  transferred  from  the  fins  to  the  air. 

All  of  the  programs  require  some  modification  to  facilitate  interfacing 
with  the  data  base  management  system  described  in  section  2.2.7. 
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2 PWB’S  8x  4 Inches 


Figure  2.  Air  Cooled  Printed  Wiring  Board 


It  became  obvious  at  the  outset  of  the  study  that  the  major  difficulty  in 
developing  an  efficient  optimization  algorithm  was  the  magnitude  and  in- 
efficiency of  the  automated  wire  routing  process.  This  process  usually 
takes  several  hours  of  computer  central  processor  time  and  then  still 
requires  some  manual  intervention  to  complete  the  circuitry.  The  wire 
routing  processes  determine  two  major  costs  which  must  enter  into  the  trade- 
offs required  for  optimum  design.  Even  one  pass  through  the  wiring  routing 
process  is  expensive.  To  make  repeated  passes  to  define  an  optimum  layout 
is  so  expensive  and  time  consuming  that  it  is  untenable  for  a production 
design  process. 

Further  complicating  the  wire  routing  problem  is  the  fact  that  the  process 
is  not  a well  defined  or  well  behaved  function  with  respect  to  reliability. 
For  example,  a slight  rearrangement  of  components  may  have  little  or  no 
effect  on  the  number  of  layers  of  PWB  or  on  the  cost  of  the  wire  routing 
process.  On  the  other  hand,  the  same  slight  rearrangement  on  a different 
board  may  have  a very  significant  effect.  As  an  example  of  this,  the  trial 
optimization  of  a PWA  such  as  is  shown  in  figure  2 indicated  that  changing 
the  layout  from  the  best  wiring  layout  to  the  optimum  reliability  layout 
increased  the  cost  by  109%  if  the  same  number  of  layers  were  retained  (see 
Table  ll.  Allowing  some  additional  flexibility  in  wire  routing  reduced  the 
cost  increase  to  90%.  For  these  two  layouts,  the  reliability  increase 
changed  from  23%  improvement  to  14%.  In  other  words,  this  reliability  gain 
was  cut  almost  in  half  with  only  a 10%  reduction  in  design  cost. 

The  above  results  might  lead  to  several  deductions.  However,  if  each 
board  is  considered  separately,  different  conclusions  become  apparent. 

For  the  first  board,  the  difference  in  cost  between  the  optimum  layout  and 
the  modified  optimum  layout  is  9%,  whereas  the  failure  rate  increased  by 
6.7%. 

The  second  board,  however,  shows  a change  in  cost  of  10.5%  with  an  increase 
in  fa  lure  rate  of  8.5%. 

Based  on  this  very  small  sample,  it  would  seem  that  there  may  be  some 
correlation  between  increasing  reliability  and  increasing  cost.  However, 
it  also  appears  that  any  derivation  from  the  best  electrical  wiring  layout 
causes  a sharp  (100%)  rise  in  cost.  The  same  sort  of  behavior  is  evident 
when  the  number  of  layers  is  allowed  to  increase,  but  there  is  a slightly 
lower  maximum  cost. 

For  purposes  of  optimization  trade-off  calculations,  the  wire  routing 
routine  should  consider  at  least  four  different  layouts.  They  would  include 

1)  Best  wiring  layout 

2)  A perturbation  of  (1)  in  the  direction  of  10%  improved  reliability 

3 ) Opt i mum  re  1 i ab i 1 i ty  1 ayout 

4)  A perturbation  of  (3)  in  the  direction  of  10%  reduced  reliability 
(improved  routeabi 1 i ty . 
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BOARD 

PART 

NUMBER 

CONSTRAINT 

SET 

f AILURE 
HATE 

NUMBE  R 

OF 

LAVERS 

DESIGN 
PRICE  S 

644040/ 

Optimized 

bi 

mm 

6100 

644R40  7 

Optimized 

3800 

644R407 

Semi  Optimized 

ES 

LH; 

6200 

64  4 M 40  ’ 

Semi  Opt imi/e«1 

12 

3800 

644M40H 

Optimized 

71 

10 

/600 

644 R 408 

Optimized 

71 

12 

4800 

644  R 408 

Semi  Optimized 

77 

10 

6800 

644R80B 

Sem«  •Optimized 

77 

12 

4800 

&44H40  7 

Not  Optimized* 

ESI 

10 

2970 

644M40H 

Not  Optimized* 

mm 

10 

38  70 

*The  Original  Depgn  was  clone  to  keep  the  Power  in  each 
Row  Equivalent 

78  0610  VA  1 


The  design  cost  of  each  of  these  layouts  for  at  least  two  different  number; 
of  PWB  layers  must  be  predicted.  It  is  completely  impractical  to  completely 
wire  route  a PWB  this  many  times.  For  this  reason,  a simplified  procedure 
that  is  able  to  predict  the  relative  costs  of  various  layouts  without 
actually  completing  the  wire  routing  process  is  needed.  Investigations 
(primarily  i' i scussions)  with  a leading  corporation  in  the  automated  wire 
routing  business  indicate  this  is  a viable  approach.  In  the  following 
discussion,  it  is  assumed  that  such  an  estimation  technique  will  be  used 
in  the  ITAD  system.  It  is,  of  course,  recognized  that  there  may  be 
several  possible  perturbed  layouts  (items  2 and  -1  above)  which  will  give  the 
same  reliability  predictions  and  have  various  impact  on  wire  routing. 
Depending  on  the  efficiency  of  the  developed  estimation  software,  a search 
can  be  made  to  attempt  to  find  a "best"  layout  for  the  same  perturbed 
re  1 i ab  i 1 i ty . 

The  LCC  optimization  block  diagram  is  shown  in  figure  3.  This  diagram 
is  shown  as  though  data  were  passed  directly  from  one  program  to  another 
although  it  is,  in  actuality,  passed  via  the  data  base  management  system 
! DRM ) discussed  in  detail  in  section  2.2.7.  An  attempt  to  represent  this 
is  shown  in  figure  4 where  the  circles  represent  data  in  the  DBM. 

The  optimization  routines  can  be  initiated  by  starting  either  at  the  existing 
computer  program  OPTEMP  which  calculates  and  optimizes  reliability  by 
inputting  the  PWA  part  number,  cooling  technique,  and  sink  temperatures. 

Other  required  information  is  obtained  from  the  data  base  management  system. 
OPTEMP  rearranges  the  component  layout  to  obtain  minimum  failure  rat e . 

During  this  process,  arrangements  are  stored  at  3%  improvement  increments 
so  that  when  optimum  placement  is  attained,  "snapshots"  of  approximately 
10%  subopt imal  layouts  can  be  retained  for  wire  routing  cost  trade-offs. 

These  layouts  are  transferred  to  the  wire  routing  cost  program  which 
calculates  the  number  of  layers  of  PWB  circuitry  required  and  associated 
costs  of  doing  the  wire  routing  and  artwork. 

The  wire  routing  cost  program  then  determines  the  best  layout  for  ease 
of  wire  routing  and  associated  costs.  All  of  this  data  is  stored  in  the 
DRM  and  the  best  wiring  layout  data  is  passed  to  OPTEMP  where  reliability 
calculations  arc  made.  Then  OPTEMP  perturbs  this  layout  until  reliability 
has  been  improved  by  approximately  10%.  This  fourth  layout  is  returned  to 
the  wire  routing  cost  program  for  costing  information. 

At  this  point  a table  has  been  established  in  the  DBM  as  shown  in  Table  2. 
This  data  allows  R & D costs  and  fabrication  costs  of  the  PWB  at  eight 
different  points  (eight  designs).  By  adding  the  costs  of  parts  and 
connectors,  the  investment  cost  is  established. 

The  ESC  for  the  four  layouts  can  be  calculated  using  the  modified  AH  C 
LSC  model.  Then  adding  the  various  costs  yields  an  eight  point  LCC  latice 
from  which  the  lowest  LCC  is  picked. 


Figure  3.  Optimization  Loop 
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A sensitivity  analysis  is  now  performed  to  determine  relative  sensitivity 
of  LCC  to  PWB  acquisition  cost,  parts  acquisition  cost,  PWB  R £ 0 cost,  or 
reliability.  At  this  point  the  only  parameter  that  can  be  further  optimized 
in  a reasonable  manner  is  the  parts  on  the  PWB.  Thus,  if  the  sensitivity 
analysis  shows  LCC  to  be  most  sensitive  to  either  part  cost  or  reliability, 
lower  quality  parts  can  be  exchanged  for  the  better  quality  parts.  If 
reliability  is  more  sensitive  than  part  cost,  then  better  quality  parts 
can  be  used. 


With  these  changes,  the  optimization  cycle  returns  all  the  layouts  with 
the  new  parts  to  OPTEMP  where  new  PKA  reliability  calculations  are  made 
and  the  table  updated.  The  wire  routing  cost  program  is  bypassed  since  the 
part  change  does  not  affect  the  wiring.  If  a completely  different  part  is 
substituted  land  this  would  be  done  manually)  which  does  change  the  routing, 
the  entire  cycle  must  be  repeated.  The  LCC  cost  matrix  and  sensitivity 
is  perfonned  and  compared  with  the  previous  minimum  LCC  to  determine  if 
an  improvement  has  been  made.  This  process  continues  until  no  improvement 
is  shown.  Then  a new  optimum  layout  is  defined  and  the  entire  cycle 
repeated . 

This  algorithm  does  not  necessarily  provide  an  absolute  minimum  LCC  since 
there  may  be  some  particular  layout  which  will  give  a slightly  better  result. 
However,  this  algorithm  will  give  a near  minimum  LCC  with  a few  iterations. 
Better  optimization  can,  of  course,  be  obtained  by  using  5%  reliability 
increments  and  more  data  points  in  the  latice  at  the  sacrifice  of  more 
computat ion. 


An  alternate  technique  which  has  the  potential  of  coming  closer  to  an 
absolute  minimum  but  may  require  more  or  less  iterations  is  to  do  a 
gradient  search  based  on  the  sensitivity  analysis  starting  first  at  the 
best  wiring  layout  until  no  improvement  is  shown.  Then  start  at  the  optimum 
reliability  layout  and  reduce  reliability  until  no  improvement  is  shown. 

Then  compare  the  two  best  LCC  designs  obtained  from  the  two  starting  points. 
Again  there  are  many  layouts  possible  to  obtain  the  same  reliability  and 
some  increment  must  be  used  so  this  still  does  not  guarantee  an 
absolute  minimum. 

The  optimization  loop  can  be  integrated  with  the  aircraft  LCC  optimization 
routine  to  make  aircraft  ECS  design  decisions.  This  can  be  done  by- 
including  the  avionics  optimization  loop  within  the  aircraft  optimization 
loop.  Thus  different  optimum  avionics  designs  and  I.CC's  which  can  be 
compared  with  the  ECS  LCC  to  find  the  best  integrated  design. 

An  approximation  of  this  technique  which  would  reduce  the  cost  but  not 
yield  as  good  an  optimum  would  be  to  not  change  the  original  PWB  layouts. 
This  would  then  only  require  recalculating  reliabilities  for  the  already 
established  layouts.  R £ D costs  and  bare  board  acquisition  cost  would 
remain  constant.  The  most  desirable  layout  may  change  and  the  part  quality- 
will  change  with  changes  in  coolant  parameters  i.e.,  flow  rate,  temperature, 
and  pressure  drop  to  provide  a new  LCC  to  compare  with  the  aircraft  LCC 
at  each  iteration. 
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It  is  obvious  that  the  optimum  avionics  design  considered  by  itself,  is 
one  where  the  aircraft  ECS  is  supplying  air  at  a high  rate  and  cold 
temperature  which  will  make  the  reliability  completely  insensitive  to 
layout.  This  would  allow  the  least  expensive  wire  routing  and  the  minimum 
number  of  layers  with  a very  low  failure  rate.  From  the  overall  weapon 
system's  point  of  view,  the  restriction  of  not  altering  the  layout  is  not 
generally  or  strictly  valid.  However,  within  a range  of  values  reasonable 
for  the  aircraft  ECS  performance  and  LCC,  the  approximation  may  be 
reasonable. 

2.2.2  Thermal  Environment  Software 
2. 2. 2.1  OPTEMP 

The  existing  computer  program  OPTEMP  which  has  been  developed  at  the 
Flight  Dynamics  Laboratory,  Wright  Patterson  AFB,  will  provide  the  means 
for  optimizing  component  positions  on  the  printed  circuit  boards.  The 
existing  version  of  OPTEMP  will  optimize  both  the  component  locations  and 
the  spacing  between  components  for  coplanar  cooled  printed  circuit  boards. 
Additional  features  which  will  be  added  to  this  program  are: 

1)  Calculation  of  coefficients  F and  & . 

2)  Calculation  of  PWB  equivalent  temperature  (this  is  defined  below) 

3)  Identification  of  components  via  circuit  symbols 

4)  Capability  of  handling  additional  cooling  schemes 

a.  In-line  conduction 

b.  Natural  convection 

c.  Forced  convection 

d.  Heat  pipes 

e.  Thermoelectric  devices 

5)  Choosing  best  cooling  technique 

6)  Orienting  PWB  based  on  equivalent  temperatures 

7)  Interface  with  DBM 

The  coefficients  F and  o<  used  in  the  expression  W = will  be  calculated 

for  each  component  based  on  two  data  points  obtained  from  RELCOMP,  i.e., 

o<  = ( 1 ) Ln  and  F = ^/^T> 

12  - 

where  2 is  evaluated  and  T2  and  /L1  is  evaluated  at  Tj . 

Tj  and  T2  are  input  by  the  user  for  RELCOMP  failure  rate  calculations. 

Components  on  each  printed  circuit  board  will  be  arranged  in  groups  of 
constant  ©< . Normally  only  2 or  3 values  of  a<  will  be  encountered.  An 
equivalent  temperature  will  then  be  determined  for  each  group  of  components 

by  TEQ  = ( >Sum/  FSu„><* i ST  slnk  - \ 
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where  A sum  = sum  of  failure  rates  of  all  components  having  ex'  o<  j 
calculated  at  a sink  temperature  . 


The  physical  meaning  of  the  equivalent  temperature  is  given  as  the  answer 
to  the  following  question.  If  all  components  were  at  the  same  temperature, 
what  temperature  would  yield  the  same  failure  rate  as  that  calculated  using 
the  actual  temperature  distribution?  That  temperature  is  the  equivalent 
temperature  calculated  as  previously  described.  In  addition  to  having  a 
physical  meaning,  the  equivalent  temperature  provides  a means  of  treating 
an  entire  printed  circuit  board  as  a component,  i.e., 


ABD  = sum£ 


EQ  + 


T T 

( sink  - 1) 


where  Tj  = sink  temperature  at  which  Tp-g  was  evaluated  and  = new  sink 

temperature . 


If  a printed  circuit  board  has  components  with  3 different  values  for  ex'  , 
then  ABD  = ( ABD)  cx  j + ( ABD)  CX.  2 + ( ABD)  cX, 

The  OPTEMP  program  currently  has  an  equivalent  temperature  routing.  However, 
it  is  based  on  the  O coefficient  of  the  first  component  read  into  the  input 
data. 


The  addition  of  circuit  symbol  numbers  for  component  identification  has 
already  been  added  to  OPTEMP  and  was  used  for  the  Task  II  effort.  This 
feature  requires  the  entry  of  the  circuit  symbol  along  with  other  component 
data  such  as  power  dissipation  and  size.  This  will  be  modified  so  that  only 
the  PIVB  part  number  need  be  entered  and  all  required  data  will  be  obtained 
from  the  DBM. 

The  ability  to  optimize  printed  circuit  board  layouts  with  cooling  schemes 
other  than  coplanar,  will  be  incorporated  into  OPTEMP.  In-lyie  cooling  is 
a common  technique  in  which  components  are  placed  on  heat  conducting  metal 
strips  which  conduct  the  component  heat  to  the  ends  of  the  board  which 
interface  with  a heat  sink.  The  temperature  profile  of  the  conductors  will 
be  calculated  by  superposition  whereby  the  contribution  of  each  component 
is  summed  with  the  baseline  conductor  temperature.  This  is  shown  graphically 
in  Figure  5 and  6 for  two  cooling  schemes.  Figure  5 is  a general  case 
of  a board  with  heat  sinks  at  each  end  which  are  at  different  temperatures. 
Figure  6 is  a simplified  case  of  only  one  heat  sink.  Optimization  of  the 
in-line  boards  will  be  achieved  by  the  modified  bubble  sort  technique 
described  herein. 

Natural  convection  cooling  will  also  be  added  to  OPTEMP.  This  cooling 
technique  is  often  preferred  because  of  its  simplicity  but  is  not  often 
practical  for  boards  with  substantial  heat  dissipation. 

It  is  assumed  that  the  p.c.  board  will  be  oriented  in  a vertical  direction 
in  an  unbounded  fluid  medium  and  that  dissipation  will  be  unifotmi. 

Components  will  be  closely  spaced  and  of  uniform  size  and  power  dissipation. 
This  allows  the  use  of  the  experimental  correlation  of  the  local  Nusselt 
number  of  a vertical  flat  plate  wHh  uniform  heat  flux  as  described  in 
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Reference  2/  and  2/.  Depending  upon  whether  the  induced  convection  currents 
are  laminar  or  turbulent,  one  of  the  following  equations  will  be  employed. 

NU„  =0.60  (Grx*  Pr)  !/5 


105  < Grx*  < 1011 
(Laminar) 


NU„  = 0.17  (Grx*  Pr)  !/4 


2 x 1013  < Grx*  Pr  < 1016 


where 


hx 


NU  = Local  Nusselt  Number  = — 

X K _ 

Grx  = Modified  Grashoff  Number  = & ^ * — 


k/2 


Pr  = Prandtl  Number  = MC  /K 

P 


h = Local  film  coefficient  (BTU/hrft2oF) 
x = Distance  from  bottom  of  board  (ft.) 


K = Thermal  conductivity  of  fluid  (BTU/hrft°F) 

g = Gravitational  constant  (ft/hr2) 

B = Volume  coeeficient  of  expansion  (1/°R) 

V = Kinematic  viscosity  of  fluid  (ft.2/hr.) 

M = Absolute  viscosity  of  fluid  (lb.  sec. /ft.2) 

Cp  = Specific  heat  of  fluid  (BTU/lb.  °F) 

The  local  film  coefficient,  component  area,  and  component  heat  dissipation 
are  all  that  is  needed  to  determine  the  component  case  temperature,  i.e., 

Tc  “ Too  - 4 • 
ha 

Component  junction  temperature  will  be  calculated  based  on  the  junction  to 
case  thermal  resistance  which  is  obtained  from  DBM  or  by  user  input  and  can 
be  obtained  from  MIL-SPEC's  or  from  vendor  catalogs.  This  technique  assumes 
that  the  board  is  in  an  unbounded  fluid  medium  and  neglects  any  turbulence 
which  might  be  generated  due  to  spacing  between  components  or  due  to 
variations  in  component  heights.  User  judgement  will  be  necessary  for  the 
case  (considered  unlikely)  of  a series  of  closely  spaced  boards. 

For  forced  convection,  the  mass  flow  is  defined  and  is  input  by  the  user. 

The  junction  to  case  thermal  resistance  obtained  from  sources  described 
for  free  convection.  The  case  to  ambient  thermal  resistance  will  be 
determined  by  = (1/ha)  (273)  (°C/W)  where  h is  the  film  co- 
efficient (BTU/hr  ft2  °F ) and  a is  the  component  area  exposed  to  the 


1/  HEAT  TRANSFER  DATA  BOOK,  General  Electric  Corp.,  Section  504.2,  pg.  14 
2/  HEAT  TRANSFER,  Holman,  J.  P.,  McGraw  Hill,  New  York,  1976. 


cooling  area  (in-).  The  constant  273  is  a unit  conversion  factor. 
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The  film  coefficient  h is  detemined  by  evaluation  of  the  dimensionless 
Nusselt  number  Nu  = (hl/k)  (1/12)  where  L component  length  measured  in 
the  direction  of  air  flow  (inches)  and  K = thermal  conductivity  of  air 
( RTU/hr  ft  °F ) . 

The  value  of  Nu  depends  upon  the  values  of  the  Prandtl  number  (Pr)  and 
the  Reynolds  number  (Re).  The  dimensionless  Prandtl  number  has  previously 
been  defined  as  Pr  - Cp/K.  The  dimensionless  Reynolds  number  is  defined 
as  Re  = (pvl/A.)  (1/12) 

where  p air  density  (LBm/ft*5) 

C - ail*  velocity  (ft/sec) 

L component  length  as  previously  defined 
u - air  viscosity  (LBm  ft  hr) 

The  equations  for  the  Nusselt  number  are  (Ref.  2' 


Nu  = 

0. 

664  pr 

1/3  Re  1/ 

Re  £ 5 

x 105 

Nu 

Pr 

1 * ( , 

,037  Re*8 

-850)  5 ? 

< 105  *- 

Re  < 

107 

Nu 

Pr 

,1/3 

1 

.227  Re 
(log  Re) 

- 2.58  I - 

. i 

850 

107  < 

Re 

Optimization  of  component  placement  is  accomplished  by  the  existing  OPTEMP 
rout i ne . 

Printed  circuit  boards  which  are  configured  with  heat  pipes  for  cooling 
do  not  require  optimization  of  component  placement  because  component 
temperatures  are  unaffected  by  position.  However,  provisions  will  be 
added  tr  OPTEMP  by  which  to  calculate  component  temperatures  and  failure 
rates  so  that  this  cooling  technique  can  be  compared  with  others. 

Component  junction  temperatures  will  be  calculated  by 

T.l  = tH/P  * q comp  (®C-H/P  + ®JC) 

. q TOTAL  a comp 

wuit  H/P  sink  * hA  condenser  * hA  evaporator 

The  user  may  input  values  for  the  hA's.  If  no  value  is  input,  a default 
value  of  1500  RTU/hr  ft-  °F, a commonly  used  "rule  of  thumb"  used  to  size 
heatpipes,  will  be  used  for  the  h of  both  condenser  and  evaporator  but  the 
areas  must  still  be  input. 

Two  cases  have  been  conceived  wherein  thermoelectric  devices  might  be  used 
and  the  capability  of  handling  these  cases  will  he  contained  in  OPTEMP. 

Case  1 considers  a printed  circuit  board  with  a cooling  scheme  resulting 
in  a few  components  which  are  too  hot.  The  basic  cooling  scheme  can  be 
maintained  and  the  problem  components  can  be  lowered  in  temperature  by  the 
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use  of  thermoelectric  devices.  Such  a device  can  suppress  the  temperature 
locally,  hut  in  the  process  it  generates  considerable  heat  which  adds  to 
the  total  heat  load  of  the  board.  These  factors  can  be  accommodated  by 
redefining  the  component  Junction  to  case  thermal  resistance  (M.l  and  the 
component  heat  dissipation  (q)  as  follows:  ( ( 0 j > (q!  T<j*g[))  q* 

q*  q ♦ q TED 

Where  the  subscript  "TED"  refers  to  properties  of  the  thermoelectric  device. 
The  value  of  the  temperature  suppresion  ( Tyj.^)  is  to  be  input  by  the  user. 
A default  value  of  20°C  will  be  included.  Similarly,  qTpp  must  In-  input, 
but  a default  value  of  1.2  x q will  be  used.  These  values  are  chosen 

because  they  are  a common  near  optimum  design  point  for  TED's.  Ref.  - • 

The  sum  of  the  two  q's  will  then  be  used  as  the  equivalent  component 
dissipation.  Optimization  of  component  locations  can  proceed  based  upon 
that  applicable  for  the  main  cooling  scheme  with  the  redefined  vnl  es 
accounting  for  the  thermoelectric  device. 

Case  2 considers  an  entire  board  configured  with  a thermoelectric  device 
(or  many  small  ones)  with  the  hot  side  of  the  thermoelectric  device  attached 
to  a constant  temperature  sink,  for  this  case,  the  component  temperature 
will  not  be  a function  of  position  so  optimization  is  not  necessary. 
Temperature  will  be  calculated  by  Tj  ( ^ 0 ) (q)  4 (80s)  (q)  ♦ 

*Tsink  “ tTK1)'  where  H jc.  is  the  component  junction  to  case  thermal 

resist  ance 

Bcs  is  the  component  case  to  sink  (thermoelectric  device)  thermal 
resistance 

q is  the  component  heat  dissipation. 

The  relationship  of  OPTEMP  with  other  programs  in  1TAD  and  the  Data  Rase 
Management  System  (DRM)  is  showin  in  figure  7.  Initial  data  on  the  cooling 
parameters  and  the  PWA  to  be  analyzed  are  input  via  punch  cards  or  inter- 
active remote  terminal.  Component  data  is  obtained  from  the  DRM.  Two 
temperatures,  T;  and  T.,  (40°C  and  90°C  unless  otherwise  specified)  art' 
used  by  RE1.C0MP  to  calculate  two  values  of  ,V  for  each  generic  type  of 
component  on  the  assembly.  These  are  used  to  determine  K and  *■*  for  the 
component  type.  OPTEMP  then  optimizes  the  PWR  layout  using  the  cooling 
scheme  designated  or  searches  for  the  best  cooling  scheme  it'  required  data 
is  provided.  The  data  associated  with  the  chosen  optimum  placement  is 
passed  to  the  DRM  for  updating  and  for  use  by  the  wire  routing  program 
(PCPRA)  or  the  routing  cost  estimating  program,  or  the  l.CC  program  if 
parts  type  optimization  is  being  done. 
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2.2 .2 .2  Optimization  Algorithms 

The  optimization  routine  used  in  OPTEMP  to  arrange  components  is  a 
straightforward  "bubble  search"  routine  amply  described  in  reference  4/  . 
This  algorithm  was  designed  specifically  for  the  coplanar  cooling 
technique  described  but  is  applicable  to  forced  and  natural  convection 
cooling  schemes.  In  all  of  these  schemes,  except  unbounded  natural 
convection,  the  air  temperature  at  the  component  being  examined  is  a 
function  of  only  the  sum  of  the  dissipations  of  upstream  components. 

The  temperature  of  the  k^  component  can  thus  be  expressed  as 
k . 

Tk  = To  >7^—  * 0 qk  + eca  \ 

i ■>  mCp  jc 

where  Tq  is  air  inlet  temperature 

0jc  is  junction  to  case  resistance  of  component  k 
0ca  is  case  to  air  resistance  of  component  k 
i = 1 is  at  the  air  inlet  end  of  the  PWB 


This  can  be  expressed  as 

_1_ 

¥ ftiCn  qi 


Tk  = To  + 9ja  Ik  - ACp^,' 

In  the  case  of  unbounded  natural  convection,  this  can  be  replaced  by 

qk 

Tk  = To  + ejc  qk  + hTxJ  A 


i.e.,  the  temperature  of  the  component  is  a function  only  of  its 
location  and  power. 

In  the  case  of  conduction  cooling,  the  present  algorithm  cannot  be  used 
because  the  temperature  of  each  component  is  affected  by  the  position  of 
every  other  component.  Referring  to  figure  5,  the  temperature  at  any 
location  is  the  sum  of  the  temperature  rises  at  this  location  caused  by 
all  the  other  components  at  other  locations  and  the  temperature  rise 
caused  by  the  component  located  at  this  position. 

The  temperature  rise  at  location  i due  to  a component  to  the  right  of 

i in  figure  5 (components  are  numbered  from  right  to  left)  can  be  expressed 

as: 

ATij  = (jjJr  Rjj  ) qi  !li  J Uj 

Rjr  + Rjl  Rjl 

1/  COMPUTER  AIDED  THERMAL  DESIGN  OF  AVIONICS  FOR  OPTIMUM  RELIABILITY  AND 
LIFE  CYCLE  COST,  Dr.  Arnold  H.  Mayer,  Wright-Patterson  AFB , March  1978. 
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R i 1 = resistance  from  i to  left  sink 

R,r  = resistance  from  i to  right  sink 

Rjl  = resistance  from  j to  left  sink 

Rj r = resistance  from  j to  right  sink 

qj  = heat  dissipation  of  component  i 

The  temperature  at  i due  to  components  to  the  left  of  i can  be  expressed  as 
_ , R jr  R j 1 s Rir 


Tij  =(Rj7-Aj1)  qj  R~7  : 1 < J 

In  the  case  where  the  left  and  right  sink  temperatures,  Tj  and  Tr,  are  not 
equal,  another  term  is  required. 

(Tr  - Tl)  Ril 

Ril  4 Ri r 

Thus,  the  temperature  at  i,  Tj , is  expressed  as: 


Rjr  Ril 


j = 1 Rjr  4 Rjl  J 


(Tr  _ Tj)  Ril 

Ti  = Ti  * 4 

1 * R • i R • 

"ll  4 ir 


£ Rji.  Rlr  q. 

j = i 4 1 Rjr  4 Rjl 


For  a single  sink  on  either  end,  this  expression  reduces  to: 

Ti  = Tl  4 RU  ±_  qj  4 £ 

J = 1 J = 1 4 1 ll 

If  component  positions  are  numbered  from  left  to  right  then 
n i - 1 


V Ti  ♦ Rn  + rjl  qj 


The  failure  rate  is  expressed  as: 

*\  n,  ‘ kTk 

A = £ FUe 

k = 1 K 

and  the  object  is  to  minimize  A . 
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Initially,  an  attempt  was  made  to  use  a dynamic  programming  approach  to 
assure  a global  optimum  solution.  This  has  been  done  and  is  explained 
in  detail  in  Appendix  B.  However,  this  approach  requires  substantially 
more  computer  time  (in  fact  for  more  than  10  components  it  is  impractical) 
and  there  is  an  intuitive  feel  that  an  appropriate  bubble  sort  routine  will 
yield  comparable  answers  even  if  the  rigor  of  dynamic  programming  is  not 
present . 

The  sorting  technique  to  be  incorporated  in  OPTEMP  for  conduction  cooled 
PWA's  involves  simply  renumbering  the  components  in  order  of  increasing 
temperature,  whatever  their  location,  as  shown  in  figure  8.  The  bubble 
sorting  technique  can  now  be  used  by  comparing  the  sum  of  the  failure  rate 
of  components  one  and  two  in  their  original  and  switched  positions.  If 
the  total  A-  is  lower  when  switched,  they  are  left  that  way.  Then  1 and  3 
are  compared  in  the  same  manner  and  the  process  continues  until  the  last 
component  is  compared  with  its  adjacent  component. 

The  process  can  continue  until  completion  in  two  ways.  A complete 
pass  can  be  made  and  then  the  components  again  renumbered  for  the  next 
pass.  This  process  continues  until  no  switches  are  made  in  the  last 
pass.  This  technique  is  illustrated  in  figure  9 and  10. 

A second  method  is  to  renumber  the  components  each  time  a switch  is  made 
and  work  backward  switching  and  comparing  all  components  until  component  1 
is  reached.  Then  the  comparison  can  resume  at  the  location  of  the  original 
switch. 

The  latter  method  will  usually  converge  faster.  However,  because  the 
objective  function  is  monotoni cal ly  increasing,  the  algorithm  will  always 
converge.  The  optimization  can  be  terminated  when  gains  become  insignifi- 
cant or  allowed  to  continue  until  all  gains  are  zero.  The  maximum  number 
of  passes  should  be  limited  to  N (N-l)  although  convergence  should  occur 
long  before  this  limit  is  reached. 

2. 2. 2. 3 SINDA 

In  addition  to  the  capabilities  to  be  included  in  OPTEMP,  the  ability 
to  model  certain  thermal  control  techniques  associated  with  avionics 
equipment  will  be  available  in  S1NDA.  The  thermal  control  techniques 
required  for  ITAD  are  electrical  heaters,  thermostatically  controlled 
flow  devices,  and  variable  conductance  heat  pipes. 

Electrical  heaters,  constant,  time  dependent,  or  temperature  dependent 
are  easily  simulated  by  several  options  within  SINDA.  The  "SIT"  and  "SIV" 
options  used  in  the  heat  source  definition  block  allow  respectively  time 
and  temperature  varying  source  interpolated  from  any  table  input  into 
the  array  data  block.  Alternately,  the  "D1DEG1"  subroutine  can  be  used 
in  the  VARIABLES  1 block  to  describe  varying  heat  sources. 

The  simulation  of  thermostats  is  alsoineluded  as  a SINDA  subroutine 
called  THRMST.  This  routine  returns  a "1"  or  a "0."  These  can  be  used 
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III.  PERFORM  ONE  PAaS  THROUGH  COMPONENTS  NUMBERED  IN  THIS 


Figure  9.  Optimization  Routine 


THERE  WERE  ONE  OR  MORE  SWITCHES,  GO  TO  H 


Figure  10.  Optimization  Flow  Diagram 
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either  as  "flags"  to  turn  on  and  off  a heater,  or  can  be  used  with  a 
multiplier  to  give  direct  calculation  of  the  heater  or  flow  controlled  b\ 
the  thermostat.  A description  of  this  subroutine  is  included  in  Appendix 
A. 


Variable  Conductance  Heat  Pipes  (VCHP)  can  be  simulated  in  SINDA  by  using 
the  temperature  dependent  conductance  option.  Using  the  "SIV"  conductor 
option,  the  condenser  area  can  be  varied  with  condensate  temperature. 

The  simulation  of  phase  change  cooled  systems  is  also  a standard  subroutine 
in  SINDA.  Actually,  three  subroutines  are  available  to  simulate  expendable 
ablation,  expendable  boiling,  and  non-expendable  melt  and  freeze  types 
of  phase  change  material.  These  subroutine  descriptions  are  also  included 
in  Appendix  A. 

Evaluation  of  temperatures  achieved  with  air  having  free  moisture  will 
be  accomplished  by  an  iterative  solution  to  the  First  Law  of  Thermodynamics . 
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|^(wcpT)2  - (wcpT) J 

j>cpT)2  - (wCpT)^ 
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q ■ heat  input  (BTU/hr)  GIVEN 

m " mass  air  flow  (Lbm/min)  GIVEN 

cp  - air  specific  heat  (BTU/lb  °F)  GIVEN 

T2  " exhaust  temp  (°F)  UNKNOWN 

■ inlet  temp  (°F)  GIVEN 

''LIQUID  " ™LIQUId/‘^AIR  UNKNOWN* 

wVAPOR  ■ ^apor/^air  unknown* 

cPLIQUID  " specific  heat  of  liquid  (BTU/lb  °F)  GIVEN 

CpVAPOR  " sPecific  heat  of  vapor  (BTU/lb  °E)  GIVEN 

hf_  ■ latent  heat  of  lporization  (BTU/lb)  GIVEN 
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P.  - (Pc)(10) 


- L 

T 


a + bx  + cxJ 


1 + dx 


where  a,  b,  c,  d,  ?c,  Tc  are  constants 

x - Tc  - T 

P6  - VAPOR  PRESSURE 

T - TEMP  (°K) 

Pv  " Ps  WLIQUID  ^ 0 
WVAPOR  = .622/ (P/Py  -1) 

WLIQl'ID  “ WU20  ‘ WVAPOR 

*CAN  RE  CALCULATED  BY  SE1  ARi\TE  EQUATIONS 


2.2.3  Dynamic  Environment  Software 

The  capabilities  of  a number  of  computer  programs  were  analyzed  relative  to 
their  applicability  to  determine  the  vibration  environments  at  equipment 
boundaries  and  the  equipment  response.  All  programs  that  were  analyzed 
are  in  the  public  domain  of  available  programs.  A summary  of  the 
capabilities  of  the  candidate  programs  along  with  other  comments  is 
given  below. 

SAP- IV  - SAP-IV  is  a general  purpose  finite  element  program  that  performs 
statis,  modal,  transient,  frequency  response  and  spectrum  analysis.  The 
program  has  an  extensive  library  of  elements  and  has  both  modal  and  prop- 
ortional damping.  "The  unique  feature  of  SAP  IV  is  that  it  is  the  only 
general  purpose  program  that  is  simple  enough  so  that  users  have  tried 
to  modify  it  for  their  needs. 5/  Additionally,  SAP-IV  is  currently  resident 
at  the  WPAFB  computer  facility  at  the  Avionics  Lab. 

SUM-10,  SHOCK  AND  VIBRATION  CONTROLS  COMPUTER  PROGRAM  REVIEWS  AND 

SUMMARIES,  Edited  by  Walter  and  Barbara  Pikley,  1975,  Shock  and 

Vibration  Information  Center,  U.  S.  Department  of  Defense. 
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One  difficulty  encountered  with  SAP-IY  is  that  it  has  no  random  vibration 
response  analysis  capability.  This  difficulty  can  be  overcome,  however, 
by  usingg^he  frequency  response  option  of  SAP-IY  with  a suitable  sub- 
routine - that  calculates  the  output  power  spectral  densities  (PSD's) 
and  RMS  response  of  specified  degrees  of  freedom  for  a given  input  PSD. 

NASTRAN  - A large  versatile  general  purpose  finite  element  program,  NASTRAN , 
contains  extensive  linear  dynamic  analysis  capabilities.  Due  to  the 
versatility  of  the  program  and  the  bulk  of  the  documentation,  considerable 
experience  is  required  for  its  use.  In  addition,  NASTRAN  is  very  in- 
efficient (costly)  for  small  problems.  NASTRAN  Level  15  is  available  at 
WPAFB  on  the  CDC*  6600. 

FLExZ  - A WPAFB  computer  program  for  vibration  prediction  of  fighter 
aircraft  equipments,  FLEX  is  currently  available  on  the  CDC  6600  and 
Cyber  73  computers.  This  program  predicts  vibration  spectra  at  equipment 
mounting  locations  on  fighter  aircraft  due  to  boundary  layer  (aerodynamic) 
inputs  for  straight  and  level  flight,  buffet  turns,  take  off,  and  landing. 

This  program  is  very  useful  in  predicting  the  vibration  environments  for 
a large  class  of  aircraft.  In  addition,  several  future  work  areas  were 
listed  which  means  the  program  will  be  expanded  and  improved  and  that  it 
should  be  continually  supported  and  maintained  by  WPAFB  personnel. 

5 / 8/ 

PURDUE  PROGRAM—*  — - This  program  estimates  the  response  of  a linearly 

elastic  system  under  random  or  deterministic  steady  state  excitation.  An 
approximate  method  based  on  relaxation  techniques  (conjugate  gradient 
method)  is  used.  Natural  frequencies  and  normal  modes  are  not  computed. 

To  use  this  program,  the  mass  damping  and  stiffness  matrices  are  cal- 
culated outside  the  program  and  input  into  it.  The  most  obvious  use  of 
the  program  is  to  fill  the  void  in  the  SAP-IY'  capabilities  due  to  its 
lack  of  a random  analysis  capability. 

9/ 

RSPC—  - A WPAFB  computer  program  for  the  random  response  analysis  of 
structures  subjected  to  multiple  random  loads.  To  use  the  program 
eigenvectors,  natural  frequencies,  generalized  forces,  masses,  and 
damping  are  computed  outside  the  program  and  then  read  into  the  program. 

-y  VIBRATION  CONTROL  FOR  AIRBORNE  OPTICAL  SYSTEMS,  AFFDL-TR-76-145 , WPAFB . 

COMPUTER  PROGRAM  FOR  VIBRATION  PREDICTION  OF  FIGHTER  AIRCRAFT  EQUIPMENT, 
AFFDL-TR-77-  , WPAFB. 

A NEW  METHOD  FOR  PREDICTING  RESPONSE  IN  COMPLEX  LINEAR  SYSTEMS  II, 

J.  L.  Bagdanoff,  et  al , Journal  of  Sound  and  Vibration,  Yol.  53,  1977, 
pp  459-469. 

-/  BUFFET  DYNAMIC  LOADS  DURING  MANEUVERS,  AFFDL-TP-72-46 , WPAFB. 
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No  user's  manual  or  other  documentation  could  be  obtained  for  this  program 
although  it  was  possible  to  get  a listing  and  program  deck.  It  was 
discovered  during  the  course  of  this  investigation,  however,  that  certain 
portions  of  the  program  had  not  been  successfully  debugged  and  that 
caution  should  be  exercised  in  using  the  program. 

SAP  V - SAP  V is  basically  SAP  IV  with  added  capabilities.  Primary  among 
these  additional  capabilities  are  the  inclusion  of  a bandwidth  minimization 
routine,  a graphics  package,  and  rigid  elements.  In  addition,  a 1978 
objective  of  the  SAP  V Users  Group  is  the  addition  of  a random  response 
analysis  option. 

The  recommended  choice  for  a general  purpose  FEA  (finite  element  analysis) 
program  is  SAP  V due  to  its  small  size,  ease  of  use,  and  good  static  and 
dynamic  capabilities.  The  proposed  inclusion  of  a random  response 
analysis  package  remedies  the  deficiencies  of  the  previous  SAP  IV 
program . 

NASTRAN  has  capabilities  that  exceed  those  seen  as  necessary  under  the 
present  requirement  of  an  optimum  design  oriented  FEA  program. 

It  is  also  recommended  that  the  FLEX  program  be  included  in  the  ITAD 
software.  The  FLEX  program,  while  being  limited  to  one  particular  class 
of  vehicles  (i.e.,  fighters)  represents  a significant  contribution  in 
the  determination  of  the  onboard  fighter  vibration  environment  to  be  used 
as  input  to  the  avionics  system. 

Rationale  for  the  Adjustment  Of  the  Applications  Environment  Factor 

The  failure  rate  model  from  MIL-HDBK  217B  ( / =TTj  TTq  Tfp  ( t 

+ C2  Il\)  contains  a temperature  dependent  factor  ff t and  an  application 
environment  factor  /T\-  The  n't  factor  is  well  defined  analytically, 
and  its  value  can  be  minimized  thru  the  optimum  placement  of  components 
for  maximum  reliability.  The  application  environment  factor,  or 
as  it  is  called,  is  a general  factor  based  on  vibration  levels  and  other 
parameters  such  as  sand,  dust,  humidity,  and  temperature  cycling.  The 
breakdown  of  what  percentage  of  the  /7“L  factor  is  due  to  each  parameter 
is  not  clearly  defined  though  some  work  in  this  area  has  been  undertakenlH/ 
and  is  presently  being  studied  at  RADC . 

As  the  current,  state  of  the  art  exists,  the  exact  contribution  of  the 
structural  reliability  to  the  overall  reliability  calculation  is  unknown, 
and  much  more  study  and  testing  is  needed  in  the  area  of  the  structural 
reliability  as  it  relates  to  the  overall  system  reliability  figure. 

It  is  obvious  that  a P.W.  board  or  other  piece  of  electronic  equipment 
that  has  been  weel  designed  structurally  to  minimize  vibration  effects 
will  be  more  reliable  than  one  that  has  been  poorly  designed. 

10/  DEVELOPMENT  OF  HIGHLY  RELIABLE  SOLDER  JOINTS  FOR  PWBs,  Westinghouse 


Final  Report  on  Contract  No.  NAS8-21233,  1968. 


For  this  reason  it  is  assumed  that:  (1)  the  structural  design  of  a P.W. 
board  affects  the  total  board  reliability  and  should  therefore  be  optimized 
as  much  as  possible  and  (2)  that  this  reliability  is  quantifiable.  That 
is,  that  it  makes  up  some  unspecified  fraction  of  the  g factor  (14% 
based  on  Grumman  study)  that  remains  to  be  determined  through  future  study. 
For  the  present,  it  is  sufficient  to  say  that  the  effect  exists,  and  to 
arbitrarily  pick  a percentage  (say  15%)  as  being  typical  and  then  proceed 
from  there. 

To  quantify  the  failure  rate  due  to  the  vibration  environment,  it  is 
first  necessary  to  identify  the  board  failure  modes  and  the  factors  that 
affect  these  failure  modes.  Some,  but  by  no  means  all,  of  the  mechanical 
failure  modes  are: 

° P.W.  board  dynamic  stresses  and  fatigue 
° Component  electrical  lead  wire  stresses  and  fatigue 
° Thermal  expansion  between  the  board  and  lead  wires 
° Solder  joint  thermal  stress  and  fatigue 
o Plated  through  holes  thermal  stress  fatigue 
Board  Failure  Modes 

The  dynamic  bending  stresses  in  a printed  wiring  board  is  determined  at 
the  resonance  by  considering  the  resonant  mode  shape  of  the  circuit  board, 
its  transmi ssibi 1 ity  and  the  input  load  condition.  From  these  factors  the 
maximum  bending  moment  in  the  board  can  be  determined  and  hence  the  dynamic 
stress.  From  this  data  and  the  board  material  properties  the  fatigue 
life  of  the  board  can  be  calculated. 

Component  Electrical  Lead  Wires  Stresses 

Electronic  component  parts  such  as  resistors,  capacitors,  and  flatpack 
I.C.’s  are  often  mounted  to  P.W.  boards  by  their  electrical  lead  wires. 
Failures  can  occur  in  these  lead  wires  due  to  the  large  deflections  and 
bending  stresses  that  are  developed  as  the  board  bends  back  and  forth. 

Every  time  the  P.W.  board  bends,  it  forces  the  lead  wires  to  bend  since 
the  wires  are  normally  soldered  to  the  board.  The  greater  the  PW  board 
deflection,  the  higher  the  bending  stresses  in  the  lead  wires. 

Stresses  in  the  component  lead  wires  can  be  determined  by  examining  the 
geometry  of  the  electronic  component,  the  fundamental  resonant  frequency 
of  the  boards  and  the  input  forcing  function.  The  relative  slope  of  the 
electric  lead  wires,  where  they  are  soldered  to  the  circuit  board,  will 
not  change  as  the  board  vibrates.  That  is,  the  wires  remain  perpendicular 
to  the  board.  The  wires  also  remain  perpendicular  to  the  component  body 
where  the  lead  joins  the  body  as  the  body  of  the  component  is  much  stiffer 
than  the  lead.  Therefore,  as  the  circuit  board  bends  back  and  forth  the 
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electrical  lead  wires  will  bend  as  shown  in  Kig.  11.  From  these  curves 
the  stress  in  the  lead  wires  and  the  solder  can  be  determined  by  analyzing 
the  lead  wires  as  a beam.!^/  * ' ■ 

Thermal  Expansion 

Due  to  the  mismatch  of  the  expansion  coefficients  between  the  P.W.  boards 
and  the  metals  used  in  the  component  leads  and  solder  thermal  induced 
stresses  arise  in  the  solder  joints.  These  stresses  can  cause  cracking 
of  the  solder  joints  and  failure  of  the  components  due  to  cold  solder 
joints.  The  stresses  due  to  thermal  cycling  of  solder  joints  has  been 
studied  extensively  and  there  exist  analytic  means  to  calculate  these 
st  resses 


Fatigue  of  Solder  Joints 

Even  when  solder  joints  are  properly  made,  high  bending  or  shear  stresses 
coupled  with  the  stress  reversals  due  to  vibration  can  lead  to  fatigue 
failure.  Many  solder  joint  fatigue  failures  are  due  to  a "shear  tear-out" 
type  of  failure  induced  by  the  bending  of  the  electrical  lead  wires.  This 
type  of  condition  as  well  as  others  can  be  analyzed  and  then  automated .i£/ 


Factors  that  Affect  the  Failure  Mechanisms 


Knowing  the  potential  failure  mechanisms  that  can  occur  helps  the  designer 
to  avoid  possibly  major  problems  that  can  occur  due  to  the  component 
layout  on  the  boards.  It  is  necessary,  however,  for  the  board  designer 
to  be  aware  of  the  factors  that  affect  the  failure  mechanisms.  A few 
of  the  primary  factors  will  be  discussed  below: 

Probably  the  greatest  factors  that  affect  the  mechanical  reliability  of 
a P.W.  board  are  the  board  orientation  and  geometry.  The  orientation  of 
the  hoard  determines  the  dynamic  loads  that  the  board  will  see.  For  this 
reason  it  is  logical  that  the  plane  of  the  board  should  be  coplanar  with 
the  two  axes  that  have  the  highest  vibration  input  loads  and  normal  to 
the  axis  with  the  smallest  dynamic  input  whenever  feasible.  The  board 
geometry  will  determine  the  fundamental  mode  shapes  of  the  board  which 
in  turn  determines  the  board  deflections  and  also  the  board  and  component 
stress. 

The  material  properties  of  the  board  along  with  any  ribs,  heat  sinks, 
or  other  mechanical  stiffeners  determine  the  stiffness  of  the  board.  These 
factors  in  conjunction  with  the  board  geometry  and  the  boundary  conditions 
determine  the  natural  frequency  of  the  board. 

The  material  properties  of  the  electronic  components  leads  as  well  as 
the  solder  used  also  affect  the  thermally  induced  vibrational  stresses. 


Ii/  VIBRATION  ANALYSIS  FOR  ELECTRONIC  EQUIPMENT,  Dave  Steinberg,  John  Wiley 
R Sons,  N.  Y.,  1973 
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The  degree  of  damping  in  a P.W.  board  affects  the  transmission  of  the  input 
forcing  function  to  the  board  and  thus  has  a significant  effect  upon  the 
dynamic  load  the  board  and  the  components  experience.  Unfortunately,  the 
transmissibility  of  a board  is  difficult  to  predict  without  physically 
testing  the  board.  Steinburg  1J_/  gives  a general  rule  of  thumb  for 
predicting  the  transmissibility  at  resonance  as 

Q = K ( fn)  -5  * K <■  2 

where  f is  the  natural  frequency  of  the  board  in  Hertz  and  K is  a constant 
based  on  the  board  type  and  mounting.  Some  of  the  factors  that  affect  the 
selection  of  K are  given  below  (see  Ref.  2 for  an  extensive  list). 

o P.W.  board  natural  frequency 

o Multiple  layer  boards 

o Heat  sink  straps 

o Stiffening  ribs 

o P.W.  board  supports 

o Board  connectors 

With  an  estimate  of  the  transmissibility  and  the  board  natural  frequencies, 
it  is  therefore  possible  to  calculate  for  a known  forcing  function  the 
environmentally  induced  stresses  in  the  board  and  the  components  using 
standard  analysis  procedures.  It  is  also  possible  to  computerize  this 
analysis,  thus  relieving  the  design  engineer  of  the  tremendous  volume 
of  calculations.  The  engineer  can  thus  be  made  aware  of  the  effect  of 
changes  to  his  board  design  or  layout  and  any  potential  problem  areas. 

Method  of  Implementation 

This  section  will  deal  with  the  steps  necessary  to  calculate  the  stress 
levels  and  margins  of  safety  of  the  board  and  components  that  make  up  a 
P.W.  board  assembly. 

The  first  step  in  analyzing  a board  design  is  the  determination  of  the 
natural  frequency  and  the  mode  shape  of  the  fundamental  board  resonance 
as  in  most  cases  the  effect  of  the  fundamental  mode  usually  dominates 
the  other  board  resonance  effects.  To  determine  the  board  natural 
frequency,  certain  information  should  be  defined.  This  information 
consists  of  the  following.  Board  Geometry  - The  board  geometry  is 
basically  the  shape  of  the  given  board  - its  length,  width  and  thickness, 
plus  any  additional  structure  such  as  stiffener  ribs,  heat  sinks,  multiple 
layers,  etc. 

Beard  Material  Properties  - The  board  material  properties  is  a broad  term 
that  describes  not  only  the  composition  of  the  board  itself'  but  also  the 
properties  of  any  ground  places,  conformal  coatings,  stiffening  ribs, 
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etc.  that  will  contribute  to  the  overall  stiffness  of  the  board  (the  effect 
of  the  stiffness  of  any  components  will  be  neglected  in  determining  the 
natural  frequency).  Needed  for  the  different  materials  that  comprise  the 
board  are  the  modules  of  elasticity  (E),  Poissons  ratio  (\f ) and  the 
material  density  (^). 

Board  Weight  - The  weight  of  the  board  must  be  known  to  determine  the 
natural  frequency.  This  weight  consists  of  the  weight  of  the  board  due 
to  the  density  of  the  material  making  up  the  board  as  well  as  the  weight 
of  any  electrical  components  or  connectors  that  are  mounted  on  the  board. 

Ir.  practice  it  will  be  assumed  that  the  total  weight  of  the  board  is 
uniformly  distributed  over  the  board  area  and  that  the  location  of 
individual  components  has  no  effect  on  the  board  natural  frequency. 

Boundary  Conditions  - The  boundary  conditions  specify  whether  an  edge  of 
the  board  is  fixed  (can't  translate  or  rotate)  free  (can  both  translate  and 
rotate)  or  simply  supported  (can  rotate  but  not  translate)  or  some  combin- 
ation of  all  three  boundary  conditions.  An  accurate  representation  of 
the  board  boundary  conditions  is  critical  in  determining  both  the  resonant 
frequency,  mode  shape  and  stress  levels  in  the  board  and  the  components. 
Therefore,  the  boundary  conditions  must  be  determined  accurately  as 
even  the  expected  input  level  can  affect  the  boundary  condition.  For 
example,  a board  edge  that  can  be  assumed  clamped  under  a low  vibrational 
input  level  may  become  simply  supported  under  a higher  input  level. 

Once  the  above  information  has  been  specified,  it  is  possible  to  calculate 
the  natural  frequency  and  mode  shapes  of  the  P.W.  board  using  one  of  two 
procedures.  The  first  procedure  is  to  model  the  board  as  ar.  assembly 
of  finite  elements  and  then  use  the  SAP  V computer  program  to  determine 
tie  mode  shape  information.  The  SAP  V program  will  determine  tie 
normalized  board  displacements  and  rotations  at  each  mode  point  ir  the 
model.  This  information  can  be  saved  for  later  use  in  deteimining  the 
board  stress  and  moments. 

The  second  method  consists  of  using  analytic  approximation?  based  on 
the  previous  given  information  and  an  assumed  mode  shape  that  satisfies 
the  boundary  conditions  to  obtain  the  natural  frequency  of  the  board. 

The  Raleigh  method  of  equating  the  board  strain  energy  to  its  kinetic 
energy  is  one  method  of  obtaining  the  board  natural  frequency.  The 
assumed  mode  shapes  used  in  the  Raleigh  method  would  be  based  on 
trigonometric  or  polynomial  series  approximation.  Appendix  C illustrates 
the  derivation  of  the  board  natural  frequency  based  on  the  Raleigh  method. 

By  using  the  analytical  approximation  approach,  it  is  possible  to  build 
up  a catalog  of  board  mode  shapes  and  natural  frequencies  based  on  common 
board  geometries  such  that  the  design  engineer  could  quickly  make  cal- 
culations. If  a new  board  configuration  came  up  that  could  be  analyzed 
by  this  means  its  approximation  could  then  be  added  to  the  catalog  or 
design  table.  This  approach  would  eliminate  the  need  for  performing  a 
time  consuming  finite  element  analysis  with  each  new  board  configuration. 
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With  the  knowledge  of  the  board  mode  shape  and  natural  frequency,  the 
response  of  the  board  to  a sine  or  random  loading  could  be  determined.  If 
the  SAP  V computer  program  was  being  used,  the  board  displacements  and 
rotations  as  well  as  the  bending  moments  in  the  board  itself  could  be 
obtained  directly  from  the  program.  If  an  analytic  approximation  had 
been  used  to  obtain  the  board  natural  frequency,  additional  calculations 
would  have  to  be  done  to  obtain  the  dynamic  loads  and  response  of  the 
P.W.  board.  These  calculation,  while  not  as  precise  as  those  obtained 
using  the  SAP  V program,  would  nonetheless  provide  adequate  information 
to  the  engineer  regarding  the  loads  in  the  board. 

Finally,  by  knowing  the  dynamic  loads  and  moment  distributions  on  a board, 
the  board  stress  levels  and  component  stress  levels  can  be  calculated. 
These  calculation  would  be  implemented  as  follows.  Each  component  on  the 
board  would  have  its  mechanical  properties  such  as  component  weight,  lead 
wire  material,  separation,  area,  and  area  moments  of  inertia  stored  along 
with  its  electrical  properties.  By  defining  the  location  and  orientation 
of  the  component  on  the  board  as  well  as  the  method  of  mounting  the 
component,  it  is  a relatively  straight  forward  procedure  to  calculate 
the  stresses  in  the  component  leads  and  solder  joints  using  the  previously 
determined  board  deflections  and  bending  moments. 


The  numerous  analytical  formulas  that 
component  stress  will  not  be  presented 


have  been  derived 
here  .12/  11/  12/ 


for  determining 


Proposed  Technique 


for  Adjusting 


rr. 


While  it  is  obvious  that  the  makeup  of  the  Tf g factor  is  not  known,  it 
is  proposed  that  there  be  more  latitude  given  to  the  designer  to  influence 
the  factor  than  the  fatalistic  approach  given  in  MIL-1M7B.  This  proposal 
is  based  on  the  previously  given  arguments  regarding  effect  of  board 
orientation  etc.  To  include  the  effect  of  the  board  design  in  the  selec- 
tion ol  the  //  j.-  factor  requires  that  a baseline  board  and  orientation 
be  selected  for  each  of  the  possible  environments  based  on  the  statistical 
analysis  of  experimental  data  relating  board  mechanical  failures  to  the 
environmental  condition.  With  this  baseline  board  defined,  the  definition 
of  which  would  include  the  board  orientation,  size  mounting  conditions, 
and  stiffness,  it  would  be  possible  to  account  for  the  effects  of  good 
design  practices  (and  poor  ones  depending  upon  the  design  contraintsl  to 
incrementally  affect  the  baseline  value  of  77"g.  With  this  latitude 
allowed,  however,  it  then  becomes  possible  to  optimize  the  design  and 
layout  of  a board  for  both  mechanical  and  thermal  reliability.  Psycho- 
logically. the  allocation  of  credit  for  sound  design  would  .set  as  an 
incentive  to  promote  it. 


In  the  case  of  boards  using  all  the  same  type  of  components,  e.g.,  of 
identical  si/e  and  construction,  the  mecht-nieal  stress  is  the  same  on  any 


12  FORMULAS  FOR  STRESS  ANP  STRAIN,  Roak , R.  , McC.raw  Hill,  4th  Edition, 
New  York,  1965. 
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of  the  components  at  a given  board  location.  Different  components  will 
have  different  stresses  at  corresponding  locations.  Assuming  that  the 
calculated  stresses  can  be  correlated  with  values  of  some  portion  of 
the  IT e factor,  the  beard  can  be  divided  into  2 dimensional  sectors  with 
corresponding  J7\  values.  For  example,  in  the  equation  for  failure  rate 


A=  771  iTq  Ifp  (C1  TT T 


the  TT g can  be  partitioned  into 
77 e =^E/S  +//E/t 

where  //  E/S  = that  portion  of  //  g affected  by  stress 

' I E/t  = that  portion  of 77 E affected  by  other  factors 


A table  of  values  is  then  generated  relating  the  x,  y location  on  the 
PWB  with  a JT E/S  value  which  can  be  used  in  the  optimization  algorithm. 


A similar  table  is  established  for  each  type  of  component  on  the  board. 


In  the  case  where  all  components  on  a board  are  identical  in  size, 
construction,  and  materials  (this  is  rare),  the  sum  of  all  the  IT p/g 
will  be  constant  no  matter  what  layout  changes  are  made.  Thus,  only 
the  /7j^optimization  need  be  considered. 


While  only  an  outline  of  a technique  is  presented  here  with  a minimum 
of  supporting  evidence,  it  is  clear  that  a large  gray  area  exists  that 
bears  extensive  investigation. 


2.2.4  Reliability  Software 


For  the  various  reliability  capabilities  required  within  ITAD,  two 
existing  software  programs  were  investigated:  The  RADC  Optimized 
Reliability  and  Component  Life  Estimator  (ORACLE)  and  the  Westinghousc 
Reliability  Computational  Program  (RF.LC0MP).  Both  programs  are  com- 
puterized MIL-HDBK-217B ’ s and  both  must  be  modified  and  recompiled  for 
inclusion  as  an  integral  part  of  the  ITAD  program.  ORACLE  is  optimized 
for  a time-share  operation,  particularly  with  regard  to  the  input  and 
output  functions,  where  the  failure  rate  of  a part  is  sought  given  certain 
modifying  conditions.  RELC0MP,  on  the  other  hand,  is  the  more  versatile 
of  the  two  programs,  being  geared  towards  procesring  parts  indentured 
to  a module,  a printed  circuit  board,  a line  replaceable  unit  (LRU),  and/ 
or  a system.  RELC0MP  is  particularly  well  structured  to  allow  trade-off 
studies  involving  such  reliability  dependent  variables  as  part  quality 
and  cooling  provisions.  Additionally,  Westinghousc  is  naturally  already 
familiar  with  the  RELC0MP  program  and  in-house  funds  are  used  in  part 
to  update  the  program  to  the  latest  version  of  MIL-HDBK-217B  and  to 
refine  RELC0MP  to  meet  new  requirements  such  as  ITAD.  For  these  reasons, 
RELC0MP  was  chosen  for  inclusion  in  the  ITAD  program.  The  following 
paragraphs  describe  RELC0MP  in  more  detail. 


RELC0MP  provides  a flexible  method  of  computing  the  reliability  parameters 
for  electronic  systems,  both  existing  and  proposed.  As  depicted  in 
figure  12,  the  program  accepts  data  from  either  parts  lists  and  application 
data  for  existing  systems  (referred  to  by  Westinghouse  as  the  EDIT  file 
whose  development  and  maintenance  are  described  in  figure  17  or  engineering 
data  prepared  for  proposed  or  developmental  systems. 

The  program  orders  the  data  provided  (either  the  EDIT  data  or  engineering 
data)  in  part  number  order  and  compares  it  to  the  reliability  dictionary 
tape  (figure  14).  The  dictionary  tape  has  the  descriptive  data  required 
to  determine  the  failure  rate  of  each  part  type.  The  program  then 
computes  the  failure  rate  ( equal  to  the  number  of  failures  per  million 

hours)  for  each  part  application  in  accordance  with  the  component  parts 
parameters  provided  by  the  dictionary.  The  appropriate  application 
adjustment  factors  (actual  or  assumed)  are  also  applied.  This  data  is 
listed  in  part  number  order  and  indicates  where  the  part  ir  used  in  the 
system . 

The  failure  rates  are  for  the  most  part  computed  by  the  methods  of 
MIL— HDBK-2 1 7B  with  all  parts  not  covered  by  the  dictionary  tape  being 
listed  so  that  manual  failure  rate  entries  can  be  made. 

The  part  failure  rates  arc  combined  in  accordance  with  the  system 
configuration  and  lists  all  failure  rates  from  the  component  level  to 
the  system  level  by  assembly.  The  program  will  take  the  assembly  failure 
rates  and  combine  them  in  accordance  with  the  system  reliability  model 
to  arrive  at  a system  reliability  prediction. 

For  the  purposes  of  analysis  and  trade-off  studies,  the  program  has  the 
capability  of  computing  for  any  combination  of  application  parameters 
(temperature  or  stress)  and  providing  matrices  to  illustrate  individual 
or  combined  effects  upon  component,  assembly,  system,  or  math  model. 

To  integrate  RELC0MP  with  the  ITAD  system,  it  is  only  necessary  to 
interface  it  with  the  new  Data  Base  Management  System  (DBM)  which  will 
be  used  in  place  of  the  old  EDIT  and  DICTT0NARY  files. 
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2.2.5  LCC  Software 

The  Westinghouse  approach  to  developing  an  1TAD  LCC  software  system  is 
based  on  maximum  utilization  of  existing  computer  programs.  This  approach 
is  described  in  the  following  paragraphs.  In  addition  to  defining  a 
suitable  approach,  Westinghouse  has  developed  : preliminary  ITAD  LCC 
model.  This  model,  which  is  based  on  the  AFAL  KAK  I.SC  model  with 
developments  by  Westinghouse,  is  also  described  and  illustrated  in  this 
section.  This  preliminary  ITAD  LCC  model  was  used  as  the  primary  tool 
to  quantify  LCC  for  the  design  trade  studies  in  task  II  of  this  program. 

2.2.5. 1 LCC  Approach 

Integrated  thermal  avionic  tTTADl  of  multi-layered  printed  circuit 
boards  (PCB’s)  is  intended  to  make  improvements  in  reliability  and 
reductions  in  life  cycle  costs  (LCC).  ITAD  is  a system  of  computer 
programs  being  developed  for  use  by  PCB  designers.  It  is  intended  to 
provide  an  effective  tool  to  optimize  practically  all  physical  parameters 
to  achieve  the  objective  of  maximum  reliability  and  minimum  LCC.  It 
results  in  a PCB  design  with  a thermal  optimum  arrangement  of  the  parts 
on  a multi-layer  PCB  to  achieve  the  above  objective.  This  section 
describes  a variety  of  techniques  and  procedures  being  investigated 
to  demonstrate  the  feasibility  of  the  ITAD  approach  for  the  Flight 
Dynamics  Laboratory  of  the  United  States  Air  Force . In  this  approach, 
particular  emphasis  is  directed  to  the  application  of  reliability  and 
LCC  computer  models  for  this  design  approach. 

This  section  includes  a flow  diagram  illustrating  the  relationship  and 
sequence  of  various  computer  models  necessary  to  optimize  reliability 
and  LCC.  These  computer  models  include  those  required  for  thermal 
calculations,  reliability  calculations,  and  rearrangement  logic.  The 
LCC  optimization  loop  is  described  and  includes  constraints  imposed 
by  wire  routing,  and  the  number  of  layers  in  a bare  board  assembly. 

Also  described  is  the  impact  of  the  three  LCC  major  cost  elements; 

(1)  research  and  design  (RKP),  (2)  investment  (including  production 
costs),  and  (5)  operation  and  support  costs. 

ITAD  optimizes  LCC  by  considering  that  the  reliability  of  a printed 
wiring  assembly  is  highly  dependent  on  the  temperature  distribution  of' 
the  components  on  the  printed  wiring  assembly . This  temperature  distribu- 
tion can  thus  be  modified  by  a rearrangement  of  components.  ITAD 
rearranges  the  components  in  a systematic  manner  to  yield  a minimum 
life  cycle  cost. 

In  order  to  use  ITAD  to  optimize  life  cycle  costs,  it  is  necessary  to 
complete  certain  prerequisites.  The  omission  of  any  prerequisite  step 
makes  it  impossible  to  produce  valid  life  cycle  costs.  These  steps  are 
listed  below  in  the  sequence  in  which  they  are  performed: 


1)  Circuit  Analysis 

;’ ) System  Design  Configuration 
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3)  Mechanical  and  Thermal  Stress  Analyses 

4)  System  Reliability 

5)  System  Costs 

6)  Maintenance  Engineering  Analysis 

7)  Optimum  Repair  Level  Analysis  and  Spares  Analysis 

8)  Life  Cycle  Costs 

The  relationship  between  these  steps  is  shown  in  Figure  15,  Relationship 
Between  Decision  Making  Models  and  Techniques.  Each  of  these  steps  is 
briefly  explained  in  the  following  paragraphs,  including  their  application 
in  the  ITAD  approach. 

System  Design  Configuration 

The  design  configuration  is  the  basic  starting  point.  Initially  the 
system  design  evolves  as  a block  diagram  for  all  functions  to  be 
performed  in  the  design.  ITAD  requires  a design  configuration  down 
to  the  level  of  a PCB  with  all  component  parts  identified.  For  this 
application,  it  is  mandatory  that  the  design  configuration  be  mechanized 
so  that  a computer  can  be  used  to  generate  a system  top-down  breakdown, 
and  an  alpha-numeric  part  number  listing  of  all  parts  in  the  design. 

System  Reliability 

The  computer  program  used  to  calculate  system  reliability , (either 
RELC0MP,  or  ORACLE  is  based  on  MIL-HDBK- 21 7B , Reliability  Prediction  of 
Electronic  Equipment  and  is  used  to  compute  the  system  reliability  for 
the  initial  design.  In  the  past  the  designs  for  each  PCB  were  accepted 
once  they  were  successfully  tested.  This  baseline  reliability  also 
locked  in  both  the  system  operational  availability  and  resultant  re- 
curring maintenance  costs.  The  ITAD  approach  further  improves  the  design 
reliability  of  each  PCB  with  a thermal  optimum  arrangement  of  the  parts 
on  a multi-layer  PCB. 

The  reliability  model  stores  a dictionary  of  unit  failure  rates  for 
all  component  parts  in  the  design.  The  initial  run  of  the  model  is  used 
to  identify  any  missing  failure  rates  not  contained  in  the  dictionary, 
buch  items  are  added  in  an  update  to  the  dictionary.  OEACLE  would  flag 
such  additions  for  validation. 

Many  equations  from  MIL-HDBK-217B  are  contained  in  the  reliability  model, 
including  the  failure  rate  equations  for  microelectronic  devices  such  as: 

Ap  ” ^L  ^C1  T + C2  ^ f) 
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Figure  15.  Relationship  Between  ITAD  LCC  Decision  Making 

Models  and  Techniques 


who  re : 


X V in  the  device  Luiluie  ULt* 


is 

ttq  is 

'ff  p is 

1 7t  18 
irK  is 
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tlio  device  learning  factor 
the  quality  factor 
the  device  pin  quantity  factor 
the  device  temperature  acceleration 
the  application  environment  factor 
circuit  complexity  factors 


factor 


The  reliability  model  makes  it  possible  to  rapidly  mechanise  the  compul- 
ation ot  failure  rate  predictions  for  all  PCH's  and  other  assemblies  in 
an  entire  electronic  subsystem.  This  design  tool  makes  it  feasible  to 
perform  any  quantity  of  reliability  calculat  ions  required  in  the  li'Al) 
design  optimization. 


System  Costs 

The  ITAl)  approach  requires  the  computation  of  costs  for  (1)  the  basic 
design  effort,  and  1'1'AD  design  iterations,  (2)  production,  and  t,3)  support 
spares.  The  costs  tor  IT AD  design  iterations  are  the  only  costs  which 
change  when  the  quantity  of  hoard  layers  Is  unchanged.  If  the  quantity 
ot  hoard  layers  Increases,  costs  for  production  and  support  spares  also 
increase.  Therefore,  it  is  desirable  to  accomplish  a thermal  optimization 
of  parts  without  increasing  the  quantity  of  hoard  layers.  iho  cost  of 
the  basic  design  effort  Includes  all  research,  design,  test,  and  evaluation 
tasks  Involved  In  the  basic  design  effort.  Production  and  support 
spans  costs  reflect  all  costs  for  labor  and  materials  involved  in  the 
production  effort. 

Maintenance  Engineering  Analysis 

The  maintenance  engineering  analysis  (MF.A)  Includes  the  calculations  and 
estimates  for  all  maintenance  tasks  required  to  support  the  avionic  sub- 
system. The  MF.A  includes  maintainability  elements  such  as  estimates  for 
lime.  In  hours,  for  access,  fault  Isolation,  replacement,  repair,  and  test 
of  a PCB.  It  also  includes  estimates  ol  repair  cycle  times,  in  months, 
at  both  intermediate  and  depot  levels  of  maintenance.  Applicable  labor 
rates,  in  dollars  per  hour,  for  the  Intermediate  and  depot  love  Is.  The 
MF.A  also  Includes  the  establishment  ot  the  maintenance  plan.  Phis  has  a 
direct  hearing  on  all  support  costs. 

If  the  HAD  design  involves  a change  in  the  quantity  ot  hoard  layers, 
then  related  changes,  such  as  mean  time  to  repair,  must  he  retlected  in 
the  MF.A  data. 


■17 


Optimum  Repair  Level  and  Spares  Analyses 


In  the  ITAD  approach  the  normal  Optimum  Repair  Level  Analysis  (ORLA) 
has  been  combined  with  a spares  analysis.  The  ORLA  is  a computer  program 
designed  to  make  comparative  logistic  support  cost  calculations.  The 
result  is  an  optimum  repair  level  decision  based  on  the  selection  of  the 
overall  minimum  LSC  for  (1)  repair  at  intermediate,  (2)  repair  at  depot, 
or  (3)  discard  at  failure.  Investigations  of  ORLA  decisions  for  PCB's 
on  many  different  avionic  systems  resulted,  in  almost  every  case,  of 
repair  at  depot.  Since  ITAD  involves  only  the  design  of  PCB's,  it  was 
decided  to  code  all  PCB's  for  repair  at  depot.  Thus,  for  the  ITAD 
approach  it  is  necessary  to  only  code  an  ORLA  decision  to  avoid  the  un- 
necessary expense  for  actually  running  the  ORLA  model  for  each  PCB. 

A spares  analysis  is  highly  desirable  for  the  purpose  of  deciding  the 
level  of  sparing  necessary  to  achieve  a required  operational  availability 
for  the  avionic  subsystem.  The  spares  analysis  involves  a plot  of  spares 
cost  versus  operational  availability.  Figure  lb  is  a sample  characteristic 
curve  of  these  two  parameters.  The  spares  analysis  provides  the  visibility 
of  costs  necessary  to  assure  the  required  quantity  of  spares  at  minimum 
cost.  If  the  baseline  operational  availability  is  adequate,  it  can  be 
Inserted  as  a constraint  for  the  ITAD  design.  The  result  is  a reduction 
in  spares  quantities,  and  costs,  caused  by  the  improved  reliability  of 
an  ITAD  design. 

Life  Cycle  Costs 

The  Logistic  Support  Cost  (LSC)  Model,  developed  by  the  U.  S.  Air  Force 
Logistics  Command,  was  selected  to  calculate  the  LSC  portion  of  life 
cycle  costs  (LCC) . This  model  was  modified  to  include  improvements 
developed  through  the  efforts  of  the  Avionics  Laboratory  at  Wright  P 
Patterson  P ir  Force  Base.  For  ITAD  applications,  additional  changes, 
incliding  development  costs  and  production  cost  equations,  were  added  to 
create  an  LCC  model  suitable  for  performing  ITAD  LCC  design  trade  studies. 
Thus,  this  ITAD  LCC  model  includes  the  three  major  cost  elements: 

(1)  research  and  design,  (2)  acquisition  (production),  and  (3)  support 
costs . 

This  LCC  model  generates  the  following  output  reports  at  the  option  of 


the  user: 

1J 

Input  data  file 

20 

System  variables 

3J 

LRU/SRU  (PCB)  input  data 

igure  16.  Sample  Characteristic  Curve  of  Spares  Acquisition  Cost  vs.  Operational  Availability 


4 J Spares  summary 
5^  Availability  summary 
Cost  summary 

7 j Cost  summary  by  percentage 

The  cost  summary  provides  a breakdown  of  LCC  for  each  cost  element  and 
for  each  PCB.  Thus,  it  is  possible  to  pinpoint  the  impact  of  the  ITAD 
design  on  LCC  for  both  the  PCB,  the  LRU  in  which  it  is  located,  and  the 
entire  avionic  system. 

2. 2. 5. 2 Preliminary  HAD  LCC  Model 

A sample  of  the  output  reports  for  the  preliminary  HAD  LCC  model  is 
shown  on  the  following  pages.  The  ITAD  LCC  run  for  the  optimum  design, 
low  flew,  12-layer  constraint,  for  all  SRU's  in  the  SC/PRG  LRU  and  the  DIG 
ASSY  LRU  was  selected  for  the  sample.  Each  report  includes  an  identify- 
ing heading  which  is  Identical  for  all  output  reports.  These  reports 
are  briefly  explained  in  the  following  paragraphs.  See  paragraph 

2. 2. 5. 3 for  the  definitions  of  variables  which  appear  in  the  LCC  output 
reports . 

System  variables,  values,  and  definitions  are  shewn  in  Figure  17.  The 
MTBF  divisor  in  this  case  is  that  which  is  required  to  convert  the  ori- 
ginal EAR  MTBF  predictions  to  the  ITAD  baseline  field  MTBF.  Variables 
normally  supplied  by  the  government  were  selected  from  the  same  variables 
in  the  AFLC  LSC  user  handbook.  Other  system  values  are  contractor 
estimates  for  the  latest  design  configuration  of  the  EAR  system.  Two 
additional  system  variables  were  included  in  the  input  data  but  do  not 
yet  appear  in  the  computer  output  report.  The  added  variables  are: 

1)  TS,  Total  Systems,  which  is  used  to  compute  the  operational  availability 
tor  one  system,  2)  FTMTTR,  Fetch  Time  and  Mean  Time  to  Repair,  is  also 
needed  to  compute  the  operational  availability. 

LRU/SRU  variables  for  the  SC/PRG  and  the  DIG  ASSY  are  shown  in  Figure  18. 
The  MFTBMA  values  are  the  field  mean  flying  time  between  maintenance 
actions  for  the  ITAD  application.  The  other  column  headings  are  identical 
to  the  LRU/SRU  variables  used  in  the  AFLC  LSC  model.  Two  other  LRU/SRU 
variables  have  been  added  but  do  not  appear  in  this  output  summary  of  the 
LRU/SRU  imput  data.  First,  RD,  research  and  design  cost,  in  dollars,  is 
the  total  design  cost,  including  any  applicable  design  costs.  These  data 
do  appear  in  the  cost  Summary,  Figure  21,  column  headed  R & D.  Second, 

UCM  unit  cost  multiplier  is  used  to  apply  a learning  curbe  for  a given 
production  quantity. 
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Figure  18.  ITAD  LCC  LRU/SfiU  Input  Data 


LRU/SKU  Spares  Summary  was  added  by  Westinghouse  to  provide  complete 
visibility  to  arrival  rates,  spares  allocations,  expected  backorders, 
and  weighted  MLDT.  Figure  19  is  a sample  of  this  summary  and  includes 
a breakdown  of  base,  depot,  and  condemnation  spares.  The  arrival  rates 
at  the  base  and  depot  are  indicated  by  the  columns  headed  STKB  and  STKD. 

The  sum  of  base  and  depot  spares  indicated  by  R & D,  and  condemnation 
quantities  are  shown  under  COND.  The  TOTAL  spares  quantity  is  the  sum  of 
(R  & D)  and  COND.  Spares  quantities  were  calculated  with  a Poisson 
subroutine,  and  an  LRU/SRU  expected  backorder  constraint  of  0.05.  The 
columns  headed  EXPECTED  BACKORDER  indicate  the  base  (XBOB)  and  depot  (XBOD) 
expected  backorders  resulting  from  the  allocation  of  base  and  depot  spares. 
Weighted  MLDT  is  the  LRU/SRU  prorated  portion  of  the  system  MLDT.  The 
sum  of  the  individual  LRU/SRU  MLDT  values  is  the  system  MLDT. 

The  Availability  Sunmary,  Fig.  20  is  a brief  listing  of  the  calculated 
variables,  MLDT  and  MFTBMA,  followed  by  the  system  inherent  availability, 
operational  availability,  and  operational  readiness.  ITAD  design 
optimizations  will  always  cause  an  increase  in  availability  when  spares 
are  computed  with  the  same  expected  backorder  constraint  as  used  for  the 
baseline  design. 


ivmuBiLiir  summary 


SVSUM  MTIR: 
system  not: 

SYSTEM  MFTBMA: 


1 .50 
. 3b(>9 
17.39 


INHERE  NT  AVAILABILITY:  .920575 
OPERATIONAL  A VA I L AP l T IT Y:  .903033 
OPERATIONAL  READINESS:  .9*8092 


FIG.  20  ITAD  LCC  Availability  Summary 


The  ITAD  LCC  Cost  Summary  (in  thousands  of  dollars),  Fig  21  is  a listing 
for  each  LRU/SRU  with  columns  for  the  eight  LSC  equations.  Total  LSC, 

R & D,  Unit  Acquisition  Cost,  System  Acquisition  Cost,  and  Total  System 
(SUBSYS)  LCC.  The  LRU-TOTAL  for  each  column  is  shewn  at  the  bottom  line, 
labled  LRU-TOTAL  in  the  first  column,  for  each  LRU.  There  are  two 
exceptions  which  will  be  changed  at  a future  date;  the  subtotals  for  LRU 
Unit  Acquisition  and  System  Acquisition  appear  on  the  top  line  for  each 
LRU.  The  LCC  TOTAL  SUBSYS  is  the  sum  of  LOGIST  TOTAL,  R + D,  and  SYSTEM 
ACQUIS.  For  example,  the  LCC  TOTAL  SUBSYS  for  the  DIG  ASSY  LUR  in  Fig.  21, 
is  10461.  (LOGIST  TOTAL),  plue  221.8  (R  + D),  plus  26845.3  (SYSTEM  ACQUIS) 
for  a total  of  37528.4  (TOTAL  SUBSYS).  In  addition,  the  ITAD  field  MTBF , 
calculated  in  the  model  from  the  SRU  and  chassis  MTBF's,  is  shown  under 
the  column  heading  MFTBMA.  The  last  column  shows  the  LRU  and  SRU  unit 
weights  in  pounds. 

The  ITAD  LSC  Cost  Summary  by  Percentage,  Fig  22,  shows  the  LRU  LSC  per- 
centage for  each  of  the  eight  LSC  equations,  and  the  LSC  total.  This 
table  is  useful  in  quickly  identifying  major  cost  drivers  in  an  electronic 
system. 
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Figure  19.  ITAL  LCC  Spares 
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Figure  22.  ITAL  LSC  Cost  summary  by  Percentage 


2 . 2 . 5 . 3 Detinitions  of  Variables  Used  in  HAD  LCC  Model 


The  following  list  include  all  the  variables  used  in  the  ITAD  LCC  model 
and  which  appear  in  the  LCC  output  reports  used  as  illustrations  in 
this  report. 

BA  Cost  of  new  and  peculiar  base  support  equipment  related  to  sub- 

system. 

BAA  Total  active  work  time  in  hours,  in  the  base  shop  per  moth 

BLR  Average  base  level  labor  rate  in  do Hare /manhours 

BMC  Average  base  material  cost  per  maintenance  action  expressed  as 

a fraction  of  the  unit  cost  of  the  LRU.  The  material  cost 
shall  include  the  mean  dollar  value  of  the  consumables 
utilized  in  repairing  the  LRU  and  any  subordiante  repairables. 

BMH  Average  manhours  expended  at  base  to  diagnose  and  repair,  or 

attempt  to  repair,  the  LRU  or  SRU 

BRCT  Average  base  repair  cycle  time  in  months 

BSTOC  Number  of  unique  parts  and  assemblies  which  are  to  be  stocked 

at  the  base  as  spares,  including  standard  parts 

CS  Cost  of  software  to  utilize  existing  automatic  test  equipment 

(ATE)  for  the  subsystem 

Cl  Initial  and  pipeline  spares  cost 

C2  Replacement  spares  cost 

C3  On-equipment  maintenance  or  removal  costs 

C4  Off-equipment  maintenance  cost 

C5  Inventory  entry  and  supply  management  cost 

C6  Support  equipment  cost 

C7  Cost  of  personnel  training  and  training  equipment 

C8  Cost  of  management  and  technical  data 

COND  Fraction  of  depot  arrivals  condemned 

DA  Cost  of  new  and  peculiar  depot  support  equipment  related  to 

subsystem 
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DAA 

DLR 

DMC 

DMH 

DRCT 

EBOB 

EBOD 

FIB 

FLA 

H 

LAC 

IH 

IMH 

I PC 
JJJ 
M 

MFTBMA 

MLDT 


Total  active  work  time  In  the  depot  per  month 

Average  depot  labor  rate  In  dollars /manhour 

Average  depot  material  cost  per  cepot  maintenance  action 
expressed  as  a fraction  of  the  unit  cost  of  the  LRU.  The 
material  cost  shall  include  the  mean  dollar  value  of  the 
consumables  utilized  In  repairing  the  LRU  and  any  subordinate 
repairables  of  the  LRU. 

Average  manhours  expended  at  depot  to  diagnose  and  repair,  or 
attemot  to  repair,  the  LRU  or  SRU. 

Average  depot  repair  cycle  time  In  months  for  the  LRU  or  SRU 
Maximum  allowable  base  back  order 
Maximum  allowable  depot  back  order 

Probability  that  LRU  can  be  fault  isolated  to  an  SRU  at  the 
base  level  when  the  SRU  fails  and  LRU  removal  results 

Cost  of  unique  flight  line  AGE  (common  and  peculiar)  to  support 
the  subsystem 

Nubmer  of  pages  of  depot  level  T.  0.  ' s and  special  repair 
instructions 

Initial  inventory  management  cost  to  introduce  a new  coded 
repairable  assembly  into  the  Air  Force  inventory. 

Cost  of  interconnection  hardware  to  utilize  existing  ATE  for 
the  subsystem 

Average  time  in  manhours  to  perform  corrective  maintenance  in 
place  for  the  LRU  or  SRU 

Initial  inventory  management  cost  for  consumable  items 
Total  number  of  pages  of  T.  0.  's  (H  & J) 

Number  of  bases  expected  to  be  utilized  in  the  deployment  of  the 
system 

Mean  flying  time  between  maintenance  actions  for  each  LRU/SRU 
Mean  logistics  down  time  in  hours  per  maintenance  action 


MRF  Average  manhours /maintenance  action  for  completing 

equipment  maintenance  records 

MRO  Average  manhours /maintenance  action  for  completing 

mantenance  records 

NPLR  Probability  that  an  LRU  will  be  fault  isolated  to  a faulty  SRU 

at  aircraft  due  to  BITE,  CITS,  or  other  detection  mehtods 

NRTS  Fraction  of  LRU/LRM  removals  which,  once  removed,  are 

expected  to  be  returned  to  depot  for  repair 

OS  Fraction  of  total  force  deployed  in  overseas  locations 

OST  The  number  of  days  to  order  and  ship  the  LRU  or  the  SRU  to 

the  depot 

PA  Number  of  new"P"  coded  repairable  assemblies  within  the  LRU 

PFFH  Expected  peak  force  flying  hours /month 

PIUP  Operational  service  live  of  the  system  in  years. 

PLR  Probability  that  an  SRU  (LRM)  failure  can  be  fault  isolated  and 

replaced  at  the  flight  line. 

PMB  Direct  productive  manhours/ man/ year  at  base  level  (includes 

both  "Touch  Time"  and  time  used  in  getting  to  and  from  Job, 
getting  tools,  AGE,  etc.) 

PMD  Direct  productive  manhours/man/year  at  depot  level  (includes 

both  "Touch  Time"  and  time  used  in  getting  to  and  from  Job, 
getting  tools,  AGE,  etc.). 

PP  Number  of  new  "P"  coded  consumables  within  the  IHU 

PSC  Packaging  and  shipping  costs  for  continental  locations  in 

dollars /pound 

PSO  Packaging  and  shipping  costs  for  overseas  locations  in 

dollars / pound 

QPA  Quantity  of  the  IRU/lRM  contained  in  the  next  highest 

assembly 


off- 

on-equipment 
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Recurring  inventory  management  cost  to  maintain  a repairable 
assembly  in  the  wholesale  system  (expressed  in  dollars/item/ 

year) 

Fraction  of  failures  for  the  IRU/LRM  which  can  be  repaired 
in  place  given  that  the  IRU/lRM  has  been  fault  isolated 

Average  time  in  manhours  to  remove  and  replace  the  faulty 
IRU/LRM  for  subsequent  report 

Recurring  inventory  management  cost  to  maintain  a repairable 
assembly  in  the  wholesale  system  (expressed  in  dollars/ 

item/year) 

Fraction  of  removals  of  LRU  or  SRU  which,  once  removed,  are 
expected  to  be  repaired  at  base 

Annual  base  supply  inventory  management  cost  in  dollars/item/ 

year 

Average  time  in  manhours  to  perform  scheduled  maintenance 
(including  preventive  maintenance,  preflight,  postflight, 
and  periodic  inspections  of  the  subsystem) 

Scheduled  maintenance  interval  in  flying  hours  for  the  sub- 
system (including  preventive  maintenance,  preflight,  post- 
flight,  and  periodic  inspections  of  the  subsystem) 

Average  manhours/maintenance  action  for  completing  supply 
transaction  records 


Expected  number  ot  arrivals  at  each  base 
Expected  number  of  arrivals  at  each  depot 


Coat  of  peculiar  training  per  man  at  base  level  (including  both 
instruction  and  training  materials ) 

Cost  of  peculiar  tracing  per  man  at  depot  level 

Cost  per  original  page  of  technical  documentation 

Cost  of  peculiar  training  equipment  required  for  the  subsystem 

Expected  total  force  flying  hours  over  life  cycle 


i 


TH 

TKH 

TRI) 

UC 

W 

XROB 

XBOD 

XTYPK 


Average  manhours/maintenance  action  for  completing  trans- 
portation (01111;; 

Annual  turnover  rate  t'or  base  personnel 
Annual  turnover  rate  for  depot  personnel 
Unit  cost  of  I.RU  or  SRU 
Weight  in  pounds  of  the  l.RU  or  SRU 
Computed  l.RU /LRU  back  orders  at  base 

Computed  LRU/I.RM  back  orders  at  depot  9d 

Number  of  different  I.RU's/SRU's  that  can  be  closely  approx- 
imated by  one  given  data  set.  Note:  XTYPK  is  useful  during 
the  deve 1 oprment  phase  of  a program;  when  the  system  config- 
uration is  fully  defined  XTYPK  would  equal  l. 


•HOC  t(t  lAl.' 


Figure  23.  Multivariable  Regression  Program 
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2.2.6  Statistical  Software 

The  1TAD  regresaior  analysis  requires  that  a mathematical  model  he 
formulated  for  failure  data  of  various  electronic  equipments.  This  model 
in  conjunction  with  the  other  analysis  programs  identified  in  the  I TAD 
study  will  attempt  to  provide  a quantitative  insight  into  the  mechanisms 
affecting  failures  and  to  predict  the  performance  of  these  electronic 
equipments  in  their  operational  and  non-operat ional  environments.  The 
data  from  this  analysis  should  provide  additional  tradeoff  data  in  the 
life  cycle  cost  analysis  for  1TAP,  by  highlighting  the  dependence  of 
failures  on  the  mix  and  levels  of  combined  environments. 

A computer  program  will  be  generated  to  interface  with  the  other  IT AO 
programs  and  will  use  the  RLSTKl  RI.FORI.  step  wise  regression  program. 

This  program  is  part  of  the  1MSL  library  which  is  available  at  the 
AFSC/ASD  Computer  Center  facility  at  WPAKR.  (See  figure  23). 

2.2.6. 1 Multivariable  Regression 

Development  of  the  Multivariable  regression  model  requires  two  tasks. 

The  first  task  is  to  determine  a set  of  causal  factors  which  have  a 
significant  impact  on  the  failure  rate.  Rxisting  failure  data  from  CKRT 
tests  or  from  equipment  in  the  field  can  be  analyzed  to  propose  suspected 
relationships  between  failures  and  mixes  and  levels  of  environmental 
factors,  and  the  Multivariable  Regression  program  can  be  used  t o detennine 
the  significance  of  these  relationships.  Once  this  is  done,  the  model 
thus  developed  can  be  used  to  predict  the  failure  data  of  some  equipment 
under  a specified  scenario  by  using  the  failure  model.  A stepwise 
multivariable  regression  program  will  he  used  foi  the  analysis. 


This  program  constructs  a mathematical  model  by  adding  one  term  it  a 
t ime  from  the  supplied  terms.  At  each  step  the  most  important  set  of 
terms  are  selected  for  the  model,  first  the  single  most  important, 
then  the  two  most  important  then  three,  etc.  The  geneial  form  of  this 
model  would  be: 

A’  - «o  -t  ♦,  (x,.x2,x3...)  4 a?  (xrx2.xv...)  4.  ... 
where 


A'  Is  the  predicted  number  of  failures 

a are  linear  coefficients  determined  by  Kl SYIF 

ate  general  functions  of  any  of  out  measurable  vai  tables 

( tempeiat ure,  humidity,  vi hi  at  ion,  etc). 

The  primary  term  in  this  formula  is  expected  to  be  the  failure  rate  model 
as  proposed  in  Mil  HDRK  ;’17R.  The  parameters  for  this  model  will  be 
generated  by  the  1TAD  programs  (such  as  SINDA  and  RFlfOMP)  based  on  the 
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1 
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given  scenario. 

These  parameters  will  be  available  to  the  Multivariable  Regression  Program 
thru  the  ITAD  Data  Base  Manager.  For  example,  suppose  we  have  some 
failure  data  for  a particular  microelectronic  board,  j for  scenario  1, 

2 for  scenario  2...  n for  scenario  n.  These  scenarios  could  be 
essentially  identical  boards  subject  to  a CERT  test  procedure  or  boards 
operational  in  the  field.  For  each  scenario  compute  the  dj  (Mil.  217R 
prediction  model)  where,  say,  Xj  would  be  the  temperature  from  SINDA, 
at  each  point  in  the  scenario,  Xo  would  be  the  c term  etc.  A model 
to  explain  the  failure  data  . 

*=  (A|,  X*,  ...  AM  would  be 

A1  - ao  + 4>j  (x1,x2,...xn) 

Once  aG  and  aj  were  determined  for  this  board,  a new  improved  prediction 
(over  MIL-HDBK-217B)  could  be  made  for  this  board  subject  to  some  new 
proposed  scenario  (i.e.  10  years  in  the  field  or  other). 

A second  example  could  be  the  board  layout  problem.  Let  the  j be  the 
sum  of  observe  failure  data  for  individual  components  in  a given  layout 
for  several  scenarios  and  «([  be  the  RELCOMP  predictions  for  these 
scenarios  and  the  other  «i ' s determined  (such  as  vibration  level  versus 
position  on  the  board).  Now  the  model  could  be  used  to  predict 
for  the  life  cycle  scenario. 

Other  terms  will  be  studied  which  are  not  specifically  treated  bv  the 
MIL  217B  term  and  will  be  included  in  the  model.  Our  reference  material 
indicates  that  some  major  sources  of  discrepency  between  projected 
reliability  and  that  realized  in  the  field  include,  failures  related 
to  rapid  and  slow  thermal  cycling  and  temperature  extremes  in  both 
operating  and  non  operating  modes,  equipment  maintenance,  and  equipment 
interfaces  or  combined  environments  such  as  altitude  - humidity  and 
humidity  with  cold  temperature.  These  functions  would  be  expressed: 

where  T Is  total  life  time  and  [T]is 
n partition  of  T into  identifiable 
conditions  of  MIL  217B 


A*  is  the  predicted  number  of  failures 

r 

for  the  specific  device  under  MIL  ?17R 
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- is  a factor  due  to  rapid  thermal  cycling  for  each  operating  and  non 

operating  condition  to  be  determined  empirically  as  a set  of  factor 

similar  to  n . Later,  a functional  relationship  involving,  rate  of 
c 

temperature  change,  amount  of  change,  accumulative  number  of  cycles 
duration  of  temperature  extremes,  etc.  may  be  found. 

<fi.j  - is  a maintenance  term  based  on  quality  of  maintenance  personnel, 
frequency  of  maintenance,  etc.  and  equipment  type. 

- a factor  due  to  interface  type  and  complexity 
<J>^  - other  factors  to  be  determined. 


We  expect  the  values  of  these  functions  to  be  a table  of  factors  for 
each  equipment  type,  based  on  the  conditions  of  the  scenario,  similar 
to  the  current  method  used  in  MIL-HDBK-217B  (see  Figure  24).  In  the 
future,  if  a mathematical  formula  can  be  determined  for  any  of  the 
factors,  it  would  be  substituted  for  the  table. 

In  addition,  the  part  itioning  of  factors  such  as  TT^  (from  MIL-HDBK-217B' 
into  effects  of  humidity  and  vibration  can  be  investigated  by  proper 
selection  and  partitioning  of  the  test  data.  For  example,  if  the  model 
for  microelectronic  components  were  rewritten: 


A*  » a + a. 
o 1 

\ \ "p  \ C1 

+ a2 

L_ 

"q 

"p 

"e 

C2 

— — 

L 



and  were  generated  by  using  data  under  controlled  humidity  variation  with 
other  parameters  held  constant,  the  a^  coefficient  would  be  compared  to 
the  model  generated  by  data  for  full  variation  of  all  parameters  in/^7. 

The  regression  analysis  will  have  the  flexibility  to  handle  additional 
terms  or  refined  terms  including  cross  product  terms,  such  as  humidity- 
altitude  or  humidity  - cold  temperature  as  more  data  becomes  available 
and  as  the  1TAD  analysis  proceeds. 

2. 2. 6. 2.  Calculation  of  Means  and  Variances  of  Environmental  Factors 

A computer  program  will  be  written  to  interface  with  the  DBM  and  arrange 
data  for  computation  of  means,  variances,  and  cross  correlation  coefficient 
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for  the  environmental  parameters  (such  as  temperature,  vibration, 
humidity,  etc)  for  a specified  scenario.  The  values  of  each  parameter 
will  be  obtained  from  the  ITAD  Data  Base  Manager  as  a vector,  V, 
of  parameters  at  various  points  along  the  scenario. 

For  each  variable  of  interest,  to  compute  the  mean,  take  samples  of  each 
vj  in  order  of  increasing  time  from  each  V,  a time  varying  vector  with 
components  vj  and  use  the  right  end  point  of  (vjj,  v j i ^ j ) , v j i + 1 as  the 
value  of  vj  in  the  interval,  and  (tj+i-t^l/T  as  the  interval  length, 
over  t,  the  total  simulation  time. 


For  the  mean  of  v.  use  a e 
1 1 


vi  (vj> 


tt+l  li 


— then 


Vj  i+1 


This  algorithm  does  not  require  equally  spaced  data  and  it  retains 
information  about  variable  extreme  values  (as  opposed  to  an  average  value 
over  the  interval). 

For  a variable  such  as  ambient  temperature,  a corresponding  temperature 
value  for  each  component  must  be  generated  from  the  mission  temperature 
profi  le . 

The  variance  of  would  be  calculated  as 


’/  (vj>  ■ ° (vj  - “j  (vo)V 


II 

r 


- >'o  <*/ 


where  aA  and  p are  defined  as  above. 

In  order  that  a finite  number  of  data  points  adequately  represent  a 
distribution  some  restriction  must  be  placed  on  the  data  points  used 
for  V.  For  example,  if  the  AGRFE  9 profile  for  altitude  were  used,  then: 
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at  least  the  7 points  shown  in  25  must  be  used.  In  addition,  if  the 
value  of  v-  changes  rapidly,  points  such  as  (a),  (b),  (c),  (d)  should 
be  included,  based  on  availability,  the  length  of  the  interval  (t^+1-t^) 
and  the  size  of  the  change  of  vj  from  point  i to  i+1. 

The  values  of  the  vj  for  each  component  as  well  as  /^-j,  and  tTj  for 

each  v-  will  be  available  through  the  Data  Base  Manager. 

J 

2.2.7  Data  Base  Management/Graphics 

The  data  base  management  (DBM)  software  within  the  ITAD  system  is 
designed  to  provide  the  following  general  features: 

• Simple  access  by  all  software  modules  within  the  ITAD  system 

• Flexible  interface  and  data  base  structure  to  allow  for  future 
changes 

• Easy  access  by  users  to  create  and  maintain  the  data  base 

• Take  maximum  advantage  of  library  files  to  minimize  user  supplied 

data 

The  block  diagram  of  the  data  base  management  system  is  shown  in  Figure  26. 
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The  user  specifies  a design  for  subsequent  HAD  analysis  by  means  of 
two  design  files  designated  as  DF1  and  DF  2 whose  contents  are  defined 
in  Figures  27  and  28. 

2. 2. 7.1  Design  File 

The  structural  design  file,  Design  File  1 (DF1),  specifies,  for  each 
assembly  in  a system,  the  specific  characteristics  of  that  assembly 
at  the  latest  point  in  the  design  cycle  as  well  as  information  to 
customize  each  usage  of  a component  part  of  that  assembly.  The  assembly 
characteristics: 

Configuration  specification 

Life  Cycle  cost  information 

Cost,  MTBF,  MTTR 

Drawing  information 

Lines,  circles,  test  characters  (i.e.,  all  graphics  necessary  to 
draw  the  assembly  excluding  the  component  parts) 

Routing  information 

Parts  area,  restricted  areas 

The  assembly  header  information  thus  contains  all  information  about  the 
assembly  excluding  the  component  parts.  The  drawing  information  specifies 
in  the  case  of  a printed  wiring  board,  the  board  outline,  the  various 
tooling  holes  and  text  (e.g.,  the  assembly  drawing  number  which  is  to 
appear  on  the  silk  screen). 

The  parts  list  information  specifies  not  only  the  parts  on  an  assembly 
but  also  the  customizing  information  relating  to  each  particular  use 
of  a component.  The  parts  lists  include: 

1)  Circuit  symbol  such  as  U61  and  symbol  location  and  orientation 

2)  Part  no. 

3)  X,  Y,  Z,  i.e.,  location  of  the  part  origin  with  respect  to  the 
assembly  origin. 

4)  Fixed  location  flag  (flag  set  if  part  location  fixed  by  designer) 

5)  0 orientation  of  the  part  with  respect  to  the  coordinate  system 
of  the  assembly. 


6)  Part  power  dissipation 

7 ) Part  temperature 
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Assembly  drawing  Number 

Configuration  Information 


Revision  List 


Cost  — LCC  information 

Failures 

Drawing  data 


Parts  list 


Item  CS,  P#,  XYZO,  T, 


Figure  27.  Physical  Characteristics  (DF1) 


Assembly  drawing  number 

Signal  name  - Socket  pin  list 
Signal  name  - Socket  pin  list 


Figure  28.  String  List  Data  (BF2) 
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8)  Part  failure  rate 

9)  Part  environmental  factors 


The  part  number  used  in  the  parts  lists  is  the  means  by  which  the  various 
library  files  may  be  accessed  to  obtain  additional,  generic  data  about 
a part  such  as: 


1)  reliability  factors 


2) 

Voltage  current  ratings/stresses 

3) 

Part  Thermal  resistance,  ©jc 

4) 

Part  size 

5) 

Part  weight 

6) 

Part  cost 

7) 

Part  thermal  capacity 

2. 2. 7. 2 

Library  Files 

All  library  files  are  accessed  by  part  number  and  are  designed  to  store 
related  types  of  information.  Thus  the  parts  drawing  file  contains 
graphical  data  which  when  materialized  for  a particular  part  on  a part- 
icular assembly  specifies  all  data  necessary  to  "draw"  that  component. 
The  reliability  file  contains  the  coefficients  and  other  parameters 
which,  when  combined  with  the  customizing  temperature  and  environmental 
factors  allow  failure  rates  to  be  calculated  and  subsequently  stored 
back  in  the  design  file  for  use  by  other  programs  in  ITAD. 

2.2.7. 3 Design  File  2 

The  second  design  file  specifies  electrical  information  for  each 
printed  wire  board.  The  file  is  accessed  by  assembly  drawing  number 
as  shown  in  Figure  29. 
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Assembly  drawing  # X 

Signal  name  i socket  pin,  socket-pin... 
Signal  name  j socket-pin... 


Assembly  drawing  # Y 


Figure  29  - DF2  Format 

The  data  in  design  file  2 for  a given  assembly  is  organized  by  signal 
name  or  net  and  lists  all  the  component  socket  pins  having  that  particular 
signal  name  where  socket  is  a circuit  symbol  designator  such  as  U39 
and  pin  refers  to  a pin  number  on  that  socket. 

In  order  to  explain  the  file  structure  further,  the  interfacing  mechanism 
to  OPTEMP,  wire  routing  (PCPRA)  and  an  interactive  graphics  system  will 
be  presented. 

2. 2. 7-4  OPTEMP  Interface 

The  component  placement  optimization  process  requires  P.C.  board  data, 
component  data,  and  cooling  data.  The  steps  listed  below  define  in 
general  terms  how  this  data  is  obtained  assuming  initiate  run. 


1)  Read  system  name 

2)  Read  PC  board  assembly  number:  PCN  (from  user) 

3)  Access  design  file  1 ( DF1 ) to  get  data  for  PCN 

PC  Board  Data 

1)  Read  board  dimensions 

2)  Read  routing  data  to  get  parts  placement  area 
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Component  Data 

1 ) Read  parts  list  for  PCN  to  obtain: 
a*  Part  number 

b.  Circuit  symbol 

c.  Part  location 

d.  Fixed  location  flag 

NOTE:  If  the  flag  is  not  set,  then  the  part  location  infor- 

mation, if  present,  implies  initial  location.  If  the 
flag  is  set,  the  location  must  be  present  and  pre- 
served by  OPTEMP. 

e.  Orientation 

f.  Failure  rate 

2)  Access  the  parts  graphic  library  for  each  part  to  obtain  component 
size  information 

Cooling  Data 

1)  Read  user  supplied  cooling  parameters  and  boundary  conditions 


The  results  of  an  OPTEMP  run  are  updated  component  location  information 
temperature  and  failure  rate  which  are  stored  in  DPI  overwriting  the 

initial  information. 


2. 2. 7. 5 PCPRA  Interface 

The  wire  routing  process  requires  PCB  data,  component  data  and  inter- 
connexion information.  The  following  steps  are  required: 
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1)  Read  system  name:  PF1 

Read  PC  board  assembly  number:  PCN 
3)  Access  design  file  DF1 


PCB  Data 

l)  Read  DF1  to  Ret  routing  area  information  including  restrictions 
such  as  "keep  out  areas" 


Component  Data 


l)  Read  parts  list  to  get: 

a)  Part  number 

b)  Circuit  symbol 

c)  Part  location 

d)  Part  orientation 

) Read  reliability  library  for  each  listed  part  number  to  get  part 


type  (5400,  54161,  etc.),  package  type  (DIP)  and  number  of  pins 
3)  Create  the  "BUCKET  Category"  which  relates  circuit  symbol  to  type 
(5400  1)31,  etc.)  for  input  to  PCI'KA 
•• ) Read  the  TY  PE  - MODEL  library  for  each  component  type  to  create  the 


•TYre-MODKL  Category"  which  relates  type  to  cotrr  >nettt  pin  pattern 
template.  (All  14  pin  dips  have  one  pin  pattern.  In  pin  dips  have 


another,  etc.) 

<:*AWIW  DATA  library  tor  ea.  It  model  type  to  create  the 


nit  the  pin  lot  it  ions,  with  respect  to  the 


6)  Read  design  file  2 for  PCN  to  get  the  Interconnect  information  in 
terms  of  signal  name,  circuit  symbol  and  pin  number  to  create  the 
"STRING  LIST  Category"  for  PCPRA. 

Thus  all  information  required  for  running  PCPRA  can  be  obtained  from 
the  design  files  and  the  generic  data  library  files. 

Graphic  Interface 

The  data  necessary  to  drive  an  interactive  graphics  systems  such  as 
STAGING  is  obtained  from  the  ITAD  design  and  library  files  by  the  following 
steps. 

1)  Read  system  name:  DF1 

2)  Read  assembly  number:  PCN 

3)  Access  DF1  for  PCN 

4)  Read  and  output  assembly  drawing  information 

5)  For  each  item  in  the  assembly  parts  lists  to  obtain 

a)  Part  # 

b)  Part  location 

c)  Part  orientation 

d)  Circuit  symbol 

6)  For  each  part  # read  the  graphic  library  file  and  perform  the 
following  steps: 

a)  Read  the  graphic  information  for  a part 

b)  Calculate  the  absolute  coordinates  for  each  graphical  primitive 
(lines,  circles,  etc.)  by  applying  the  translation  and  rotation 
specified  in  the  assembly  parts  list. 

c)  Output  the  "materialized"  data. 

I 
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2. 2. 7. 6 Data  Bast'  Management  Software 

A variety  of  approaches  to  providing  the  software  capable  of  handling 
the  wide  range  of  1TAD  requirements  have  been  considered.  The  study 
concluded  that  the  required  DBM  software  is  commercially  available  with 
sufficient  capability,  documentation,  availability,  and  maturity.  One 
such  system,  system  S2000  from  Management  Resources,  Inc.,  examined  in 
depth,  can  elegantly  handle  the  storage,  maintenance,  retrievel,  and 
interfacing  requirements.  SUOOO  currently  requires  5CKio  words  of  core 
storage,  but  MRI  is  currently  developing  an  overlay  version  which  should 
substantially  reduce  this  number.  The  entire  1TAD  system  will  require 
approximately  60K10  words  if  the  overlay  version  is  not  used.  If  it 
is  used,  SINDA  is  the  program  controlling  memory  size.  The  memory 
required  for  SINDA  depends  on  the  size  of  problem  for  which  SINDA  is 
used.  For  most  requirements,  40Kio  is  adequate. 

The  VENUS  data  base  management  system  currently  installed  at  WPAFB  was 
also  examined  but  found  to  be  unacceptable  since  no  provision  exists  for 
software  interfacing  with  FORTRAN  programs. 

The  second  viable  possibility  concluded  from  the  study  is  that  the  DBM 
requirements  can  be  provided  by  writing  a new  software  system  specifically 
designed  for  HAD.  In  so  doing  advantage  can  be  taken  of  known  structural 
relationships  and  interfacing  requirements,  thus  resulting  in  a more 
efficient  solution  both  in  terms  of  compactness  of  storage  and  computer 
run  time.  The  disadvantage  is  the  longer  term  maintenance  considerations 
given  the  anticipated  evolution  of  ITAD  with  associated  modifications  to 
the  DBM  and  the  significantly  higher  initial  cost.  The  current  price 
of  Sf 000  is  $45,000  without  any  discount  which  will  be  available  through 
1979.  It  is  estimated  that  the  custom  solution  will  be  five  times 
more  expensive  than  the  S11000  solution. 

Since  the  $2000  system  can  adequately  meet  the  requirements  of  ITAD 
and  has  the  flexibility  to  meet  future  requirements,  at  much  lei.:,  cost 
than  a custom  system,  it  is  recommended  that  SC  000  (available  for  the 
CDC  hardware)  be  used  as  the  ITAD  Data  Base  Management  system. 

O.C  .8  Task  I -Summary  and  Conclusions 

During  Phase  1 of  Task  I a study  was  conducted  to  determine  which  com- 
puter site  was  most  applicable  to  the  ITAD  system.  The  site  chosen  is 
the  ASD  CDC-6600  at  WPAFB.  Before  the  wire  routing  program,  PCI'RA,  at 
AFA1.  is  converted  to  the  CDC-6600,  a study  will  be  conducted  to  gain 
experience  with  the  sensitivities  of  those  parameters  affecting  ICC. 

The  parameters  affecting  ICC  which  must  be  optimized  are,  quality  of 
component  parts,  number  of  PWB  layers,  cost  of  wire  routing,  and 
reliability.  In  the  initial  (Task  II)  study  of  a coplanar  cooled  l'KA, 
based  on  708  system  with  a ten  year  life,  there  are  indications  that  the 
number  of  layers  must  he  minimized  even  at  large  increases  in  the  cost 


of  wire  routing.  This  may  not  be  true  for  systems  used  in  small  quantiti 
but  for  these  systems  LSC  gains  through  increased  reliability  will  also 
be  small. 

Whether  the  number  of  layers  is  a reasonable  optimization  parameter  can 
only  be  determined  by  trial  optimization  of  a larger  sample  of  systems 
than  has  been  undertaken  in  this  study.  An  algorithm  has  been  presented 
which  will  optimize  all  of  the  above  parameters  based  on  the  development 
of  an  economical  (in  real  time  and  computer  time)  method  of  estimating 
the  number  of  layers  and  routing  difficulty.  A larger  scale  study  might 
indicate  whether  or  not  the  number  of  layers  need  be  considered.  If  not, 
savings  in  development  of  the  wire  routing  cost  estimating  routine  will 
be  seen.  Also  a savings  in  computer  time  to  optimize  the  PWA  each  and 
every  time  TTAD  is  used  could  be  significant. 

Several  computer  programs  have  been  recommended  to  accomplish  various 
functions  of  the  ITAD  system  and  are  described  in  detail.  These  programs 
are  listed  in  Table  3.  In  addition,  where  existing  software  is  not 
available  and  existing  knowledge  and  techniques  allow,  algorithms  have 
been  defined  to  accomplish  these  tasks.  The  algorithms  presented  allow 
the  calculations  shown  in  figure  30. 

Algorithms  to  rigorously  define  wire  routing  cost  estimation  and  TT . E 
variation,  have  not  been  presented  due  to  a large  unknown  area.  The 
estimation  of  the  number  of  layers  of  printed  wiring  required  to  route  a 
PWB  and  the  calculation  of  actual  costs  to  do  this  for  various  numbers  of 
layers  has  not  been  attempted  in  a completely  automated  manner.  Although 
it  is  not  beyond  the  state  of  the  art  to  do  this,  it  was  beyond  the 
scope  of  this  study  to  attempt  to  define  this  algorithm. 

An  outline  of  a scheme  for  affecting  the  TT g factor  used  in  present 
reliability  calculations  has  been  presented.  However,  due  to  the 
uncertainty  surrounding  the  exact  nature  of  this  parameter,  a precise 
optimization  algorithm  would  be  conjecture. 
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Table  3.  RKCOMMENDED  COMPUTER  PROGRAMS 

Thermal 

S1NDA 

0PTEMP 

Structural 
SAP  V 

Reliability 

RELC0KP 

Statistical 

RLSTEP/RLF0RL 

Data  Base  Management 
MRI , S2000 

Life  Cycle  Cost 
AFLC  LSC 


Wj^re  Rout 
PCPRA 

Graphics 

STAGING 
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1.  LCC  Optimization  of  FWA's 

2.  Reliability  Optimization  of  Conduction  Cooled  FWA's 

3.  Calculation  of  Temperatures  for  IVA's  Cooled  by: 

conduction 
natural  convection 
forced  convection 
heat  pipes 

thermoelectric  devices 

4.  Calculating  Equivalent  Temperatures 

5.  Optimizing  FWA’s  with  Respect  to  Each  Other 

6.  Taking  into  Account  Free  Moisture  in  Cooling  Air 

7.  Calculation  of  Means  and  Variances  of  Environmental  Factors 

8.  Interactive  Graphics  Interface 


Figure  30.  Defined  Algorithms 


. 


80 


SECTION  111 


TASK  II-  OPTIMAL  CIRCUIT  BOARD  DESIGN  AND  COST  STUDY 


A design  and  cost  study  was  undertaken  for  a pair  of  multilayer  coplanar 
air-cooled  circuit  boards  used  in  12  locations  in  the  Electronically  Agile 
Radar  (EAR)  built  by  Westinghouse.  This  study  has  given  a preliminary 
measure  of  the  life  cycle  cost  payoff  resulting  from  the  optimization  of 
component  arrangement  on  the  printed  circuit  boards. 

3.1  THERMAL/RELIABILITY  CONSIDERATIONS 

The  construction  of  the  printed  circuit  board  pair  used  as  the  basis  for 
this  study  is  shown  in  Figure  31.  Each  board  has  six  rows  of  components. 

Up  to  sixteen  flat  pack  components  having  14  or  16  pins  each  can  be  located 
in  each  row.  Flat  packs  with  24  pins  are  somewhat  larger  so  that  rows  with 
theses  components  will  accommodate  less  than  eixteen  packages.  In  addition, 
only  four  of  the  six  rows  can  accommodate  the  larger  size  of  the  24  pin 
flat  packs. 

For  the  purposes  of  this  study,  the  inlet  air  temperature  was  held  constant 
at  85°F.  Two  different  air  flow  rates  were  considered  - 4.75  lb/min/kw 
which  results  in  an  exhaust  temperature  of  135°F  and  3.2  lb/min/kw  which 
results  in  an  exhaust  temperature  of  160°F. 

For  effective  utilization  of  the  cooling  air,  it  is  necessary  to  group  the 
components  such  that  the  power  dissipation  per  row  is  a constant  for  each 
board.  This  was  done  for  both  the  EAR  configuration  and  for  the  ITAD 
configuration.  No  attempt  had  been  made  to  specify  component  locations  for 
the  EAR  configuration.  However,  the  final  arrangement  was  fairly  good 
from  a reliability  viewpoint.  For  example,  the  408  board  which  has  18  MOS 
devices  with  a failure  rate  two  orders  of  magnitude  higher  than  most  other 
devices  was  arranged  with  these  MOS  devices  located  at  the  air  inlet  region 
of  the  board.  This  is  the  same  configuration  which  resulted  from  the  ITAD 
optimization  of  the  board.  Accordingly,  the  reliability  improvements  and 
life  cycle  cost  improvements  that  have  been  demonstrated  using  ITAD  would 
be  even  more  dramatic  when  optimizing  a poorly  configured  board. 

Components  located  on  each  of  the  two  boards  are  identified  in  Tables  4 
and  5.  Included  in  these  tables  is  An  identification  number  which  can  be 
used  to  identify  components  on  the  board  layouts  presented  later. 
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The  arrangement  of  components  for  maximum  reliability  (minimum  failure 
rate)  was  accomplished  by  the  computer  program  OPTEMP.  This  program  was 
developed  at  the  Air  Force  Flight  Dynamics  Laboratory,  Wright-Patterson 
AFB  and  furnished  to  Westinghouse  for  this  study.  Component  failure  rates 
are  calculated  by  the  equation  X- f «2l  ' which  for  the  temperature  range  of 
this  study  is  a good  approximation  of  the  temperature  sensitive  part  of  the 
failure  rate  defined  in  MIL-HDBK-217B . The  coefficients  F and  are 
constant  for  each  component.  Inspection  of  Tables  4 and  5 reveals  that 
the  value  of  F is  different  for  most  all  components  and  that  only  3 values 
of  r<  were  encountered.  One  of  the  3 values  is  associated  only  with  a 
resistor  pack  used  once  on  each  board.  Figure  32  provides  a measure  of  the 
accuracy  of  the  failure  rate  approximation  \-F  4''  for  monolithic  micro- 
electronic devices.  The  semi-logarithmic  graph  used  in  this  figure  will 
result  in  a straight  line  for  an  exponential  such  as  F 1 . The  two 

solid  curves  are  based  on  MIL-HDBK-217B . The  dashed  lines  represent  the 
approximation  used  in  OPTEMP  with  the  values  of  F and  ^ evaluated  based 
on  data  at  40°C  and  90°C.  The  approximation  is  within  15%  of  the  MIL-HDBK- 
217B  values  over  the  temperature  range  of  30°C  to  100°C.  The  majority  of 
component  temperatures  were  found  to  be  in  this  range. 

Several  component  placement  techniques  were  investigated.  These  included: 

1)  unconstrained  component  arrangement,  2)  100%  specified  component 
arrangement  for  optimum  reliability,  3)  zone  specified  component  arrangement 
for  semi-optimum  reliability,  and  4)  constant  junction  temperature  arrange- 
ment. The  last  technique  (constant  junction  temperature  arrangement)  is 
one  which  arranges  components  in  a manner  which  minimizes  the  temperature 
spread  in  each  row.  This  technique  Was  found  to  yield  a poorer  reliability 
than  the  original  arrangement.  Accordingtly , no  life  cycle  cost  studies 
were  undertaken  for  this  placement  technique. 

Complete  life  cycle  cost  studies  were  undertaken  for  the  first  three 
placement  techniques.  The  unconstrained  component  arrangement  forms  the 
baseline  for  the  life  cycle  cost  studies.  This  configuration  shown  in 
Figures  33  and  34  is  the  one  actually  used  for  the  EAR  radar  which  had  no 
thermal  placement  requirements.  The  100%  specified  component  arrangement 
represents  a mazimized  reliability  configuration  as  determined  by  OPTEMP. 
This  configuration  is  shown  in  Figures  35  and  36.  The  zone  specified 
component  arrangement  techniques  was  conceived  in  response  to  anticipated 
wire  routing  problems  with  the  100%  specified  case.  The  zone  arrangement 
technique  divided  the  board  into  zones  and  allowed  components  in  a zone  to 
be  swapped  with  any  other  component  in  that  zone.  This  was  allowed  on  a 
row  by  row  basis  only.  Component  swapping  from  one  row  to  another  was  not 
allowed  because  this  could  cause  a heat  imbalance  in  the  rows.  The  con- 
figuration resulting  from  this  technique  is  shown  in  Figures  37  and  38. 
Figures  39  and  40  illustrate  the  actual  swapping  within  the  specified 
zones . 
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Table  4.  407  BOARD  COMPONENT  INFORMATION 


r 


i». 


Identification 

Number 

Quantity 

Generic 

Number 

Power 

(watts) 

F 

od 

@ 149°F* 
(Failures  Per 
Million  Hrs . ) 

1 

2 

54S08 

.125 

.00031 

.0234 

.010 

2 

1 

54S86 

.250 

.00031 

.0234 

.010 

3 

1 

54S74 

.150 

.00066 

.0234 

.022 

4 

1 

54S04 

.114 

.00041 

.0234 

.013 

5 

2 

54S00 

.075 

.00012 

.0234 

.004 

6 

1 

54S20 

.038 

.00021 

.0234 

.007 

7 

4 

54LS04 

.012 

.00041 

.0234 

.013 

8 

1 

54S37 

.165 

.00031 

.0234 

.010 

9 

1 

54S153 

.225 

.00080 

.0234 

.026 

10 

12 

25S09 

.375 

.00143 

.0234 

.046 

11 

3 

25S10 

.525 

.00119 

.0234 

.049 

12 

3 

54157 

.150 

.00036 

.0234 

.012 

13 

1 

54LS283 

.103 

.00138 

.0234 

.045 

14 

4 

54S174 

.450 

.00138 

.0234 

.045 

15 

1 

54S151 

.225 

.00083 

.0234 

.027 

16 

1 

RP1K 

.187 

.22533 

.0045 

.441 

17 

1 

54S138 

.245 

.00076 

.0234 

.025 

18 

3 

54LS174 

.080 

.00138 

.0234 

.045 

19 

2 

54368 

.295 

.00050 

.0234 

.016 

20 

6 

54S181 

.600 

.00200 

.0234 

.065 

21 

15 

25S05 

.600 

.00250 

.0234 

.082 

temperature  related  failures  only. 


84 


Number 


Generic 

Number 

Power 
(watts ) 

F 

54LS04 

.012 

.00041 

54S74 

.150 

.00065 

55363 

.237 

.00031 

54S00 

.075 

.00012 

54S08 

.125 

,00031 

54S37 

.165 

.00031 

54S04 

.114 

.00041 

54S174 

.450 

.00138 

54LS153 

.031 

.00080 

82S11 

.600 

.01227 

54S157 

.180 

.00076 

54LS157 

.049 

.00089 

54S153 

.225 

.00080 

54LS174 

,080 

.00138 

RP1K 

.187 

.22533 

25LS09 

.055 

.00143 

25S09 

.375 

.00143 

54367 

.325 

.00050 

4060 

.500 

.00888 

f)  <a  149°F* 
(Failures  Per 
Million  Hrs . ) 


’’'Temperature  related  failures  only. 
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Figure  32.  Monolithic  Microelectronic  Device  Failure  Rate 
With  .Junction  Temperature 
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